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We calculate the light hadron spectrum in full QCD using two plus one flavor Asqtad sea quarks 
and domain wall valence quarks. Meson and baryon masses are calculated on a lattice of spatial 
size L m 2.5 fm, and a lattice spacing of a ~ 0.124 fm, for pion masses as light as m-n ~ 300 MeV, 
and compared with the results by the MILC collaboration with Asqtad valence quarks at the same 
lattice spacing. Two- and three-flavor chiral extrapolations of the baryon masses are performed 
using both continuum and mixed-action heavy baryon chiral perturbation theory. Both the three- 
flavor and two-flavor functional forms describe our lattice results, although the low-energy constants 
from the next-to-leading order SU(3) fits are inconsistent with their phenomenological values. Next- 
to-next-to- leading order SU(2) continuum formulae provide a good fit to the data and yield and 
extrapolated nucleon mass consistent with experiment, but the convergence pattern indicates that 
even our lightest pion mass may be at the upper end of the chiral regime. Surprisingly, our nucleon 
masses are essentially lineaer in over our full range of pion masses, and we show this feature 
is common to all recent dynamical calculations of the nucleon mass. The origin of this linearity is 
not presently understood, and lighter pion masses and increased control of systematic errors will be 
needed to resolve this puzzling behavior. 

PACS numbers: 11.15.Ha,12.38.Gc,12.38.Lg,14.40.-n 
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I. INTRODUCTION 

A precise calculation of the hadron spectrum is an important achievement of lattice QCD, providing a comparison 
of the experimentally known states to a first principles prediction from QCD. These calculations additionally provide 
insight into presently unexplored states, help to guide experiment, test predictions from QCD-inspircd models of 
hadrons, and set a foundation for more challenging lattice calculations. A detailed map of the hadron spectrum, 
including low lying excited states, is a challenging task for lattice QCD which is now well underway. 

In this work, we focus on the lowest lying hadrons for a given set of quantum numbers. This is part of an extensive 
program by the LHP Collaboration to study hadron structure using a mixed action that exploits the lattice chiral 
symmetry provided by domain wall valence quarks and ensembles of computationally economical improved staggered 
sea quark configurations generated by the MILC Collaboration. Since the MILC Collaboration has performed an 
extensive investigation of the spectrum using Asqtad valence quarks [1, 2], this study provides an important test 
of the systcmatics of the mixed action formulation, as well as a comparison with the spectrum determined from 
staggered fcrmions. In addition, because of the technical difficulties in the construction of baryon interpolating 
operators using the staggered quarks, this work also provides a more comprehensive study of the baryon spectrum 
than has been performed using staggered quarks alone. Since the range of pion masses treated in this work extends 
down to m T ~ 300 MeV, our lattice calculations provide a fruitful opportunity to explore the applicability of chiral 
perturbation theory for the extrapolation of hadron masses in this pion mass regime. 

In addition to this work, there now exists a large number of lattice calculations that utilize this specific mixed 
action framework. There is other work by the LHP Collaboration exploring hadron structure [3-7] , as well as related 
studies of the electromagnetic structure of hadrons [8-11]. The NPLQCD Collaboration has published a number of 
calculations focussed on hadron interactions on the lattice [12-23]. The first calculation of the hyperon axial matrix 
elements of the octet baryons was recently performed [24] and there is also a group in the process of calculating the 
kaon bag parameter [25]. Simultaneously with the development of these mixed action calculations, extensive effort 
produced the mixed action effective field theory (EFT) [26-38] , an extension of chiral perturbation theory %PT [39- 
42], which includes the relevant lattice spacing artifacts for these mixed action lattice calculations. This allows for a 
simultaneous extrapolation in the quark mass and the lattice spacing. By performing the chiral extrapolations of the 
calculated hadron masses, we will be able to determine the numerical values of the low energy constants (LECs) that 
appear as coefficients of the operators in the EFT, most importantly the physical LECs that contribute to hadronic 
processes in the continuum limit. The extracted values can then be compared with phcnomcnological determinations, 
where they exist, as well as provide predictions for those LECs which are notoriously difficult to determine without 
the ability to vary the quark masses. The real predictive power then comes when one has confidence in the numerical 
values of these universal coefficients. Once determined from one set of hadronic observables, they can then be used to 
make predictions about other observables, allowing a precision comparison with our experimental knowledge as well 
as predictions of experimentally unexplored or inaccessible physical processes. 

The rest of this paper is organized in four additional sections. In Section II we present the details of our computation 
including the tuning of the valence quark masses as well as the determination the residual chiral symmetry breaking 
inherent in finite fifth-dimensional domain-wall fcrmion lattice actions. The hadron spectrum is presented in Section III 
beginning with the meson masses and decay constants. Section III E describes the lattice calculation of the baryon 
spectrum, including the octet baryons, the nucleon, A, E and 5, as well as the decuplet baryons, the A, E*, S* and 
the In Section IV, we present the details of our chiral extrapolations. A noteworthy result of this analysis is the 
finding of unexpected and unexplained non-analytic quark mass dependence in the baryon spectrum. The behavior is 
most striking in the resulting nucleon mass calculations which are well described by a linear function in the pion mass. 
Comparing our numerical results to those of other lattice collaborations, wc find this unexpected chiral non-analytic 
behavior is also present in their lattice-data sets, suggesting this behavior is not simply the result of lattice artifacts. 
Additionally, we describe and apply the continuum and mixed action versions of heavy baryon chiral perturbation 
theory, for both two and three flavors, to our lattice calculations. All the chiral extrapolation formulae used have 
been determined previously in the literature and we refer the interested reader to the references within the text for a 
complete understanding of the derivation of these formulae. We present our conclusions in Section V. 
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TABLE I: The values of the domain-wall (DW) quark masses, the resulting pion mass obtained using the mixed action (MA) DW 
valence quarks on the asqtad MILC ensembles, together with the corresponding lightest Goldstone-boson pion mass obtained 
by the MILC collaboration to which the MA pion masses were tuned. 





0.007 


0.010 


0.020 


asqtad 

am 

u/a 

0.030 


0.040 


0.050 


# Props 

„„DWF 
am u/d 

m,K (MILC) 


3085 
0.0081 
0.1847(7) 
0.1888(2) 


3729 
0.0138 
0.2242(5) 
0.2245(2) 


3389 
0.0313 
0.3124(4) 
0.3113(2) 


2255 
0.0478 
0.3761(5) 
0.3779(2) 


324 
0.0644 
0.4338(12) 
0.4351(03) 


425 
0.081 
0.4774(10) 
0.4842(02) 



II. COMPUTATIONAL DETAILS 

We employ the gauge configurations generated by the MILC collaboration, using 2 + 1 flavors of improved staggered 
quarks using the so-call Asqtad action [43, 44]. 1 For this study we use the 20 3 x 64 lattices generated at six values of 
the light-quark masses, with the heaviest corresponding to the case of three degenerate quark masses approximately 
at the strange-quark mass and the lightest corresponding to a pion mass of m w ~ 290 MeV. The scale is set by the 
lattice spacing a = 0.12406 fm determined from heavy quark spectroscopy[2] with an uncertainty of 2%, yielding a 
spatial volume V = (2.5 fm) 3 . 

For the valence quarks, we adopt domain- wall fermions (DWF)[51, 52]; these have the arbitrarily precise realization 
of chiral symmetry required for our study of nucleon structure, and admit a straightforward interpretation of baryon 
quantum numbers, in particular parity. In order to reduce the density of small eigenvalues, the lattices are first 
HYP-smeared [53]. For computational economy, the majority of the valence propagators are computed on lattices of 
temporal extent 32 rather than 64, with Dirichlct boundary conditions imposed in the temporal direction. We refer 
to these as chopped propagators and to the those which are computed with the full length in time as unchopped. We 
improve statistics by making several computations of the spectrum for each configuration, using different temporal 
ranges. 

The physical quark fields, q(x,t), reside on the 4-dimensional boundaries with fifth coordinate 1 and L s and have 
bare quark parameter (am)® WF . The wall height M5 = 1.7 was chosen on the basis of spectral flow analyses to 
optimize the evaluation of domain wall propagators, and L s = 16 was chosen[3] to ensure that the residual mass 
characterizing residual chiral symmetry breaking is always less than 20% of the physical quark mass as discussed 
later. 



A. Tuning of the valence-quark mass 

We set the valence strange-quark mass am r/WF using the Nf = 3 ensemble corresponding to three degenerate 
sea quarks with a"i^/^ tad by imposing that the valence pseudoscalar mass equal the mass of the Goldstonc boson 

constructed using staggered quarks within the calculational precision of a few percent. This yields am^ WF = 0.081. 
We maintain the same am r/WF on the Nf = 2 + 1 ensembles. Similarly, the valence quark mass is set by tuning the 
DWF pion mass to the taste-5 staggered Goldstone-boson pion to within the calculational precision of a few percent. 
See Table I for more details. There is a freedom in the tuning choice, and for any choice, the mixed action will violate 
unitarity at finite lattice spacing. The role of the mixed action EFT analysis is to recover the proper continuum QCD 
behavior for any tuning choice [29], and our particular tuning choice provides for the most chiral valence fermions, 
which recover the QCD point in the continuum limit. For further details, see Ref. [3]. Throughout this work, we will 
refer to a lattice ensemble generated with a light Asqtad quark mass of am^ q J ad = 0.007 as the m007 ensemble, and 
similarly for the other light quark masses. 



Our results depend upon the validity of the rooting procedure used in the generation of the staggered ensembles. For recent discussion 
on this topic in the literature, we refer the reader to Refs. [45-50]. For the rest of this work, we assume this method to be valid. 
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B. Determination of m le a 

Domain wall fermions in the infinite L s limit possess an exact chiral symmetry for vanishing m q . The corresponding 
symmetry transformation is 

-> e ir5(s)e(a:) *(x,s) (1) 
*(x,a) -> §(2:, s)e- ir5 W e ^ (2) 

where ^(s) = sign( a 2 ~ 1 — s). However, at finite L s this chiral symmetry is explicitly broken by the coupling of 
left handed and right handed modes in the middle of the 5th dimension. As a result one can construct the following 
partially conserved axial vector current 

L a -1 

A^x) =- ^ r 5 (s)j M (x,s), (3) 

s=0 

where 

Ui x > s ) = \ [*( x + A, + iM+^ix, a) - s)(l - ip)U Xll &{x + A, s)] , (4) 

is the four dimensional conserved vector current that corresponds to the 4D Wilson fermion action. This current 
satisfies a Ward-Takahashi identity which in the flavor non-singlet case takes the form [52]: 

^(Al{x)0(y)) = 2m q (q(x)T a l5 q(x)0(y)) + 2(q mp (x)T a l5 q mp (x)0(y)) + i(6 a O(y)) , (5) 

where 

q{x) = ^ x ,L s -l)±^ + *(x,0)±^ . (6) 
are the physical 4D quark degrees of freedom localized at the boundaries of the 5th dimension and 

Qm P (x) = | - 1)1^ + *(*, ^)i±2* (7) 

are four dimensional fields constructed at the midpoint of the 5th dimension. The Ward-Takahashi identity of 
Eq. (5) is the same as the continuum counterpart, except for an additional term 2{q mp {x)T a ^q mp {x)0{y)) . This term 
is there only at finite L s 2 and it is a measure of the explicit chiral symmetry breaking. At long distances this term 
is proportional to 2(q(x)T a -f^q(x)0(y)) . Using the pseudo-scalar density as a probe operator 0(y) the residual mass 
is defined as 

^ t max -to ^ (q(t)T a l5 q(t)q(0)T a l5 q{0)} ' 1 ' 

where [to,t max ] is the time interval where only the ground state pion contributes to the two correlators in the ratio. 
The ratio 

R(f) = (q m p(t)T a J 5 q mp {t)q{0)T a -f 5 q(0)) 
1 ' (q(t)T a J 5 q(t)q(0)T a J 5 q(0)) ' 1 ' 

from which the residual mass is determined, is plotted in the left panel of Figure 1 for the am^' 1 " 1 = 0.01 ensemble; 
m ICS is obtained from a fit of the data to a constant. The right-hand panel shows the quark-mass dependence of m res . 
The statistical errors on the residual mass are determined using the jackknife method. 



2 For the flavor singlet current this term survives the infinite L s limit and gives rise to the axial anomaly [54, 55]. 
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FIG. 1: The left-hand figure shows the residual mass determined from the ratio R(t) in Eq. (9), for the data at am^ tad = 0.010; 
the quoted value of m re s is obtained from a constant fit to the data. Note that the effect of the Dirichlet boundary at t — 22 
relative to the source is apparent as far out as t = 14. The right-hand figure shows the quark-mass dependence of the residual 
mass. The main reason for the increase in m rcs with decreasing am^ tad is the increased roughening of the gauge field; in 
typical quenched calculations, a smaller dependence of the residual mass on the quark mass is observed. As a result, chiral 
symmetry is satisfied to a lesser degree at fixed L s as the pion mass decreases. 



III. RESULTS 



A. Fitting methodology 

An additional complication in performing fits to correlation functions constructed using the DWF valence-quark 
action arises from the oscillating terms in the DWF transfer matrix. As explained in Rcf. [56] , for the choice of domain 
wall height M5 = 1.7, lattice artifacts at the cutoff scale produce negative eigenvalues in the transfer matrix that 
cause temporal oscillations in correlation functions at short times. Since these modes at the cutoff scale decay rapidly 
and do not contribute to physical low mass states, one can fit the correlation functions used for our spectroscopic 
measurements by the usual sum of decaying exponentials plus additional oscillating exponentials e~ Mt (— 1)*. Hence, 
in this work, we have used the following form [56, 57] that includes the physical state of interest, the next excited 
physical state, and one oscillating artifact state: 

C(t) = Aoe-"^*-*^ + A 1 e~ M ^ t - t ^ + A osc (-l)*e- M <-(*-*^. (10) 

Results for the spectrum of both mesons and baryons are obtained from fits to a single correlator smeared both 
at the source and at the sink. Our fitting method is as follows. We perform fits using the form Eq. (10) introduced 
above, which enables the use of data close to the source, but the need for six fitting parameters. We also perform 
single-exponential, two-parameter fits to correlators smeared both at the source and at the sink, with a final fit range 
chosen so as to ensure a acceptable quality of fit. In general the masses are insensitive to the largest time extent 
used in the fit. We check the consistency between the two fitting methods and except where noted, obtain most final 
results using the single-exponential method. 

In general, the quoted errors are obtained using single-elimination jack-knife. However, for the study of the chiral 
behavior of the baryon spectrum, we adopt a different procedure. We chose a method that aids in accounting for the 
correlations among the fit parameters for the chiral extrapolations of the baryon masses. In the case of a calculation in 
the quenched approximation, the same number of configurations would be employed at each value of the valence quark 
mass, and therefore this could be accomplished using the jack-knife procedure above. For our full QCD calculation, 
we have different numbers of configurations at each quark mass, and therefore we adopt a bootstrap procedure in 
which we use 10,000 bootstrap samplings on each ensemble, irrespective of the number of configurations. In this case, 
the quoted errors are obtained using the 68% confidence level in the bootstrap distribution. 
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FIG. 2: The pion (left) and kaon (right) effective masses at each value of the light-quark mass, together with the single- 
exponential fits to the correlators, as described in the text. The oscillatory terms in the transfer matrix are evident in the 
effective mass close to the source. 



TABLE II: Meson masses in lattice units, calculated from the correlation functions created with the interpolating operators, 
i/jTip. Masses we were unable to extract are denoted "N/A". 



Particle 




r 


m007 


mOlO 


m020 


m030 


m040 


m050 


TV 


o-+ 


75 


0.1842(7) 


0.2238(5) 


0.3113(4) 


0.3752(5) 


0.4324(12) 


0.4767(10) 


K 


o-+ 


75 


0.3682(5) 


0.3791(5) 


0.4058(4) 


0.4311(5) 


0.4578(12) 


0.4767(10) 


ss 


o-+ 


75 


0.4827(4) 


0.4846(4) 


0.4816(4) 


0.4805(5) 


0.4818(12) 


0.4767(10) 


p 


1~ 


7» 


0.554(8) 


0.574(5) 


0.601(2) 


0.629(2) 


0.653(4) 


0.678(3) 


K* 


r~- 


7< 


0.615(2) 


0.6287(19) 


0.6397(14) 


0.6542(16) 


0.666(4) 


0.678(3) 


4> 


1" 


7< 


0.6701(11) 


0.6777(10) 


0.6774(10) 


0.6794(14) 


0.679(4) 


0.678(3) 


do 


0++ 


1 


N/A 


N/A 


N/A 


1.08(6) 


0.93(9) 


0.88(4) 


K* Q 


0++ 


1 


N/A 


1.05(6) 


0.93(3) 


0.95(3) 


0.91(6) 


0.88(3) 


/o 


0++ 


1 


0.809(14) 


0.833(13) 


0.879(18) 


0.93(3) 


0.88(7) 


0.86(6) 


Oi 


1++ 


7^75 


0.78(5) 


0.84(4) 


0.93(2) 


0.920(19) 


0.94(4) 


0.97(3) 


Ki 


1++ 


7*75 


0.89(2) 


0.915(17) 


0.960(14) 


0.954(15) 


0.96(4) 


0.97(3) 


h 


1++ 


7^75 


0.929(8) 


0.956(9) 


0.983(9) 


0.982(12) 


0.97(3) 


0.97(3) 


In 


1+- 


7i7j 


1.1(4) 


1.04(11) 


0.91(3) 


0.94(3) 


1.03(6) 


1.03(4) 


hi 


1+- 


7<7j 


0.92(2) 


0.98(2) 


0.97(3) 


0.99(4) 


1.09(13) 


1.05(10) 



B. Meson spectrum 

We use meson interpolating fields of the form 4>Tip, which overlap with the physical states listed in Table II; charge 
conjugation C applies only to particles with zero net flavor. Our single-exponential fits using the procedure described 
above are also summarized in Table II. In Figure 2, we show the effective masses of the pion and kaon, together with 
the corresponding fits. In the case of the flavor-neutral states, for example ss, <fi, etc., we only compute the connected 
contributions to the correlation functions. In Figure 3, we plot our resulting kaon, ss and vector meson masses along 
with those of the coarse (a ~ 0.124 fm) MILC calculations [1, 2]. The differences in these correlation functions are the 
direct result of differences in the valence fermion actions used in the calculations, and thus a measure of discretization 
effects. 

In Figure 3(a), one observes that the kaon and ss masses are the same within a few percent over the entire range of 
light quark masses used. This is the direct result of the tuning we have employed, discussed in Sec. II A and displayed 
in Table I. The strange quark mass was tuned on the m005 ensemble and held fixed both in all the coarse MILC 
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FIG. 3: Comparison between pseudo-scalar and vector masses with domain wall valence quarks (LHP) and staggered valence 
quarks (MILC) computed on the coarse (a « 0.125 fm) MILC ensembles. 



ensembles as well as in all the valence domain-wall propagators calculated for this work. As the pion, kaon and ss 
mesons are pseudo-scalars, the good chiral properties of the domain- wall propagators and the taste-5 staggered pseudo- 
scalar mesons protects these masses from additive mass renormalization resulting in the good numerical agreement 
observed. All other hadron masses are expected to suffer from additive lattice spacing mass corrections (beginning at 
0(a 2 )), the size of which depends on the specifics of the discretization method. Evidence of this is given, for example 
by calculations of the </> mass, which we compare to those of MILC in Fig. 3(6). However, while we observe the mass 
splitting of the is largely independent of m ff , as expected, this is not the case for the p and K* vector meson masses. 
Our calculational results tend to converge with those of the coarse MILC calculations, depicted in Fig. 3(c) and (d). 

The comparison of these meson mass calculations is complicated by the fact that they are unstable particles, at 
least at the lightest two mass points. 3 While we can not explain in detail the observed diminishing of the mass 
splitting of the p and K* masses at lighter pion masses, we do believe we understand the origin of the noted different 
behavior. There are two important issues related to the comparison of our calculation with that of MILC, which 
arise from the lack of unitarity present in both our mixed action calculation and that of MILC. These are in fact the 



3 Strictly speaking, in a Euclidean finite volume, there are not unstable particles, but rather a mixing of single and multiple particle 
states with forbidden energy level crossings [58]. Then, a particle which in an infinite Minkowski volume is unstable, in Euclidean 
finite volume may experience power law corrections to its energy levels instead of the standard corrections which scale exponentially 
as exp (— m n L). In this work we will adopt the perhaps misleading, but intuitively guided language of describing these states with 
Minkowski-space terminology. Because of the constrained kinematics in the finite volumes we work with, these would-be unstable 
particles are kinematically forbidden from decaying. Or, in more appropriate terminology, the mixing of the cigenstates of the hamiltonian 
experiencing power law volume dependence with the states experiencing exponentially small volume dependence is suppressed for the 
typical hadron correlation functions calculated on the lattices we use in this work. 
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FIG. 4: Two types of "bubble" contributions to the scalar meson: the left one is B\ in Eq. (11) while the right panel is B2 in 
Eq. (11). 

same issues which must be addressed for the calculation of the do correlator [59, 60] (which we address in more detail 
shortly). For both discretization methods, there are hairpin interactions in the decay channel, similar to Fig. 4(6) for 
do- However, in the case of staggered valence propagators, there are two additional hairpin interactions than for the 
mixed action calculation. Calculated in the appropriate EFT, these hairpin interactions are known provide a negative 
contribution to the correlation functions and become particularly relevant when the intermediate two-particle states 
are close to going "on-shell." In the continuum limit, the hairpin interactions in the two lattice calculations would be 
the same, and would vanish on the QCD line of degenerate valence and sea quarks. The second issue is related to 
the density of two-particle states near the decay threshold in the two different lattice calculations, those states which 
contribute to the p and K* correlators similar to those for the ao depicted in Fig. 4(a). In the case of our mixed-ation 
calculation, the only two-particle states which couple are those of QCD, irir for the p and Kit for the K* , however 
their masses are heavier as they are mixed valence-sea mesons [36]. In the case of the MILC calculation, there are 16 
intermediate two-particle states for each one that exists in QCD, due to the extra staggered taste degrees of freedom, 
with masses given by the known taste-splittings on the MILC ensembles [61]. It is plausible that these two issues give 
rise to the strong mass dependence observed in these vector meson mass splittings. One would expect similar issues 
with the <f>, however due to the strange quark being too heavy on these configurations, the dominant decay mode of 
the cf), the KK state, is significantly above threshold. It would have been nice to resolve this issue with a comparison 
of the fine (a w 0.09 fm) and super-fine (a « 0.06 fm) MILC calculations. However, the super-fine results are not 
available yet and the fine have a reduced strange quark mass as well as smaller lattice spacing. 

Returning to the mass of the scalar meson, we note that it has been notoriously difficult to determine, not only 
on the lattice but also experimentally. In the case of our mixed action, as mentioned above, the scalar meson is 
sensitive to taste-breaking terms in the staggered sea sector. Such an effect can be described by effective field theory 
and subtracted from the measured correlators to extract the scalar meson mass correctly in the point-source and 
point-sink case using the formulation by Prclovsck [59] . This has been performed in detail in Rcf . [62] for the case of 
domain wall valence fermions and the staggered sea fermions. There are two major "bubble" types of contributions 
(see Fig. 4) that one should include in order to extract the scalar mass (in the notation of Refs. [59, 62]): 



where 



Bi(t)= £ 



B(t) = f^(B 1 (t) + B 2 (t)). 



(11) 



k=-L 



2u4 v (m 2 



' M,i) (^w - rfi) 



(12) 



with k running over the allowed triplet integers of lattice momenta and 



E 



k=-L 



3e 



-2u) vs t 



4^ 



3e-W 
2<, 



(13) 
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FIG. 5: The left panel of this figure shows the "effective" point-point scalar meson correlator plot (as defined in Eq. (15)), 
along with the mixed-action bubble contribution from Eq. (11). Note that the symbols are the same as in Fig. 2. The right 
hand panel shows the fits to data from the lightest three ensembles where the largest "bubble" contributions dominate. 



In these equations, m vv is the pion mass calculated from the DWF valence sector, 



mui = y/A~i + 2fi stag mi , 



'T 1 



ml v + 2/j, stag mi 



l ui 



l si 



(14) 



Similarly, msi and m vs are obtained by replacing the light quark mass, mi, with the strange one. m s , in the staggered 
sector. The quantities /i s tagj Mdw, A m j x and Aj can be obtained from Rcfs. [2, 36, 61]. However, in our work, we 
only compute the Gaussian smeared propagator with Gaussian and point sink. A naive combination of these two 
correlators provides an "effective" point-point correlator, 



Cp P = Ccp{t) 2 /Cgg{P)- 



(15) 



Unlike the close match between the scalar correlator and the bubble term in Ref. [62], when compared to our "ef- 
fective" point-point correlators, the bubble-term contribution from Eq. (11) does not fully describe the unphysical 
contributions, see Fig. 5. This is not completely surprising, since our correlators use Gaussian smearing. Even though 
we project onto zero momentum in our calculation, there may be non-trivial contributions to these bubble terms. 
Therefore, we consider an alternative fit form including B(t) = diBi(t) + d2Bi(t) in Eq. (10). Unfortunately, even 
though we are able to fit all the correlators with this new formula, the large number of parameters in the fit results 
in fitted masses with large error bars. 

Similarly, abnormal behavior is observed in the strange partner Kq on the lightest sea ensemble. The Kq decays 
primarily to Ktt and only on the heaviest ensemble is it above this threshold. On the other hand, /o is the most 
stable of all the scalar mesons; there is no negative dipping behavior, and decay into two p mesons cannot happen 
here. As with the p-wave mesons, we found that the b\ masses are within the statistical error bar of the a\ meson 
on the heaviest 4 ensembles but are heavier on the lightest two ensembles. This is in contradiction with experimental 
expectations. Their strange partners, h\ and /i, are almost consistent within their statistical error bars. However, it 
is difficult for us to extract similarly accurate masses from these correlators. 



C. Decay constants 



Recently, the utility of performing chiral extrapolations in terms of m ff // ff has been demonstrated in a number of 
papers [6, 7, 12, 16, 17, 19, 20]. In addition to alleviating scale setting issues, for meson physics, the benefit of this 
type of extrapolation is understood by using the mixed action EFT in terms of renormalization and the use of "lattice- 
physical" parameters, which eliminate the leading and much of the sub-leading lattice spacing corrections [34]. Less 
restrictive but similar arguments also apply to other hadronic quantities including baryon quantities [38]. To utilize 
this method, we determine with our set of propagators following NPLQCD [16]. In addition, we use the values of 
/k/Itt in Rcf. [16] to determine fx on the m007-m030 lattices so that we can explore baryon mass extrapolations in 
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TABLE III: Pseudoscalar meson decay constants. An asterisk denotes a value of fx determined from /„- and the values of 
fa/ fn calculated in Ref. [16]. 



ID 


m007 


mOlO 


m020 


m030 


m040 


m050 




1.983(12) 


2.325(9) 


3.035(8) 


3.489(9) 


3.82(2) 


4.101(20) 




3.964(18) 


3.938(13) 


3.957(9) 


4.008(10) 


4.05(2) 


4.101(20) 


u 


0.0929(4) 


0.0963(3) 


0.1026(2) 


0.1076(3) 


0.1131(7) 


0.1162(6) 


Ik 


0.1079(4)* 


0.1087(3)* 


0.1103(2)* 


0.1122(3)* 


0.1155(7) 


0.1162(6) 



terms of m^/fx, providing insight into SU(3) breaking effects. We also calculate fx on the m040 lattice ensemble 
since this was not done in Ref. [16]. The resulting values of the decay constants are collected in Table III. 



D. Chiral extrapolation of and f„ 

One of the ultimate goals of the combined efforts of lattice QCD calculations and effective field theory is to determine 
the numerical values of the LECs of the chiral Lagrangians. This would allow us to test the range of applicability of 
chiral perturbation theory and to make predictions of obscrvables not readily computable with lattice QCD. Perhaps 
two of the simplest quantities one could imagine performing this chiral extrapolation analysis for are the pion mass 
and decay constant. Including the physical value of the pion mass, further allows for a determination of the light quark 
mass, m. The relevant two-flavor extrapolation formulae for our mixed action calculation are known to NLO [29] and 
given by, 



2Bm i 1 



2Bm (2Bm 



2Bm 



A 2 



(47T/) 2 



1 + ln 



(2Bm\ 




(16) 
(17) 



In these expressions, 1% and l\ are the renormalized LECs [40]. The unphysical coefficients Z™ and l* a arise from 
discretization effects and cannot be determined with one lattice spacing. However, at this order they act simply to 
renormalize the LO values of the chiral condensate and decay constant, B and /. The hairpin partial quenching 
parameter is simply given by (on our specific mixed action calculation) 



A 2 PQ = a 2 A l . 



(18) 



where a 2 A\ is the known coarse staggered taste splitting [61]. The mixed valence-sea meson masses are given by 



a 2 AMi: 



(19) 



where the mixed mass splitting has been calculated for the domain-wall valence fermions on the coarse MILC ensem- 
bles [36]. Unfortunately, our numerical calculation is insufficient to reliably extract the values of these LECs. For 
example, using the lightest three mass points, corrected for NLO finite volume shifts, one finds a value I3 = 4.00(5) 
where I3 is standardly defined [40], 



h = -4(4tt) 2 ^) - In 



■K,phys 



(20) 



We have not made any attempt at determining a systematic error bar to emphasize our lack of confidence in this 
number. While this number is consistent with other lattice determinations [63], in order to have confidence in the 
determination of I3 and I4 from this mixed action method, one will need calculations on at least one additional lattice 
spacing and likely lighter pion masses. 



E. Baryon spectrum 

The extraction of the baryon spectrum from a lattice calculation is obscured by the small number of double- valued 
irreducible representations admitted by the cubic symmetry of the lattice, namely G lg / U (2), H g / U {4) and G2 S / U (2), 





(«) 



FIG. 6: The left- and right-hand figures show the nucleon effective mass at am as ^ tad = 0.05 and flm^* ali = 0.007, respectively, 
from the smeared-smeared correlators. The circles denote the lattice data for the effective mass, and the solid line and dashed- 
dot line denote the single-exponential, and oscillating fit, respectively; note that lattice data is offset for clarity. 



where the g and u denote positive and negative parity, respectively, and the number in brackets denotes the dimension 
of the representation. We follow the technique introduced in Ref. [64] to construct all the possible baryon interpolating 
operators that can be formed from local or quasi- local u/d— and s-quark fields; the number of such interpolating 
fields is listed in Table IV; the G 2g / U irreducible representations are not accessible with only local or quasilocal 
fields. To yield the largest possible statistical ensemble, the correlators are averaged over the rows of the irreducible 
representation. In this paper, we present the calculations of the masses of the lowest-lying positive-parity states in 
both the octet and decuplet sectors; for the case of three flavors of degenerate quarks, SU(3) symmetry reduces the 
spectrum to single spin-1/2 and spin-3/2 ground states. 



1. Comparison of fitting forms 

We begin by a comparison of the two fitting procedures described in Section III A: "oscillating fits" of the form 
of Eq. (10), and single-exponential fits with errors derived from the bootstrap distribution. In Figures 6 and 7 
we show the effective masses for the nucleon and delta, respectively, at the SU(3)-symmetric quark masses, and at 
am u/d tad = 0-007. Also shown are the effective masses obtained from the fits using our two chosen procedures. There 
is a notable consistency between the two procedures. Whereas we will employ the single-exponential procedure in 
most of the remainder of the discussion, we show the masses obtained using fit form Eq. (10) and use the oscillating 
fit in analyzing the volume dependence of the delta mass in Fig. 14. 



TABLE IV: The number of linearly independent interpolating operators that can be formed from u/d and s quark quasi-local 
fields in each irreducible representation of the cubic group. The number in brackets denotes the dimension of the representation. 

Flavor G lg/u (2) H g/u (4) G 2g/u {2) 

N 3 1 

A 1 2 - 

A 4 1 - 

E 4 3 - 

H 4 3 - 

Q 1 2 
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FIG. 7: The left- and right-hand figures show the delta effective mass at clttl^J^ = 0.05 and am^J ali = 0.007, respectively, 
from the smeared-smeared correlators. The circles denote the lattice data for the effective mass, and the solid line and dashed- 
dot line denote the single-exponential, and oscillating fit, respectively; note that lattice data (the circles) are slightly offset for 
clarity. 



TABLE V: The upper and lower tables show ground-state masses in the octet and decuplet sectors, respectively, obtained from 
fits of form Eq. (10), as described in the text. 

asqtad 

am 

Channel 0.007 0.010 0.020 0.030 0.040 0.050 



N 



0.696(7) 0.726(5) 0.810(5) 0.878(5) 0.941(6) 0.991(5) 

0.837(4) 0.850(4) 0.886(4) 0.925(5) 0.963(6) 0.991(5) 

0.891(2) 0.902(2) 0.924(3) 0.951(4) 0.976(6) 0.991(5) 

0.784(4) 0.804(3) 0.861(4) 0.911(4) 0.957(6) 0.991(5) 



Channel 



A 

S 

n 



asqtad 
u j a 



0.007 



0.010 



0.020 



0.030 



0.040 



0.050 



0.966(13) 0.974(9) 1.005(12) 1.056(11) 1.076(14) 1.123(8) 

0.998(9) 1.026(6) 1.043(9) 1.080(10) 1.090(13) 1.123(8) 

1.053(5) 1.073(4) 1.081(7) 1.104(8) 1.104(12) 1.123(8) 

1.103(3) 1.117(3) 1.117(5) 1.127(8) 1.115(11) 1.121(9) 



2. Single- exponential fits 

Our final computations of the masses are obtained from the single-exponential fits to correlators smeared at both 
the source and at the sink, using for each flavor sector the operator that yields the smallest statistical uncertainty 
in the correlator. In order to illustrate the quality of the data on each ensemble, we show in Figures 8 and 9 the 
effective mass plots in the nucleon and delta channels for each of our u/d masses. The quality of the data degrades 
appreciably with decreasing light-quark mass, despite the increasing ensemble size, most notably in the case of the 
delta mass. One minor exception to the consistency between "oscillating fits" and "single exponential fits" occurs for 
the case of the delta effective mass at nm^' tld = 0.01, where, as seen in Fig. 9, the "single exponential fit" is visibly 
higher than the "oscillating fit" . However, as seen in Tables V and VII, the masses agree within errors, and the case 
of the delta mass is addressed in more detail in connection with Fig. 14. 

Our results for the ground-state masses for the octet and decuplet baryons are listed in Tables VI and VII, respec- 
tively. The fitting ranges are chosen so as to provide an acceptable x 2 /dof and quality of fit. The central values 
for the quoted masses are obtained from a fit to the complete ensemble. The error bars are obtained from the 68% 
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FIG. 8: The figure shows the micleon effective mass obtained from a single correlator smeared at both source and sink at each 
of the values of the light-quark masses used in the calculation; the lines show single-exponential fits to the correlator, with the 
dashed lines bootstrap errors corresponding to the 68% confidence levels on the bootstrap distributions. 




FIG. 9: The figure shows the A effective mass obtained from a singie correlator smeared at both source and sink at each of 
the values of the light-quark masses used in the calculation. The lines show single-exponential fits to the correlator, with the 
dashed lines bootstrap errors corresponding to the 68% confidence levels on the bootstrap distributions. 
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confidence level in the bootstrap distributions of the masses, using 10,000 bootstrap sub-ensembles. 

TABLE VI: The table lists the masses of the lowest-lying, positive-parity spin-1/2 octet baryons at each value of the light (u/d) 
quark mass, together with the time range used in the fit. The quoted errors are obtained from the 68% confidence level using 
10,000 bootstrap sub-samples. 



Channel 


asqtad 
am u/d 

0.007 0.010 0.020 0.030 0.040 0.050 


N 
A 
E 


16-26 0.691+ 5 

- 5 

16-26 0.780+ 3 

- 3 

16-26 0.831+ 4 

- 4 

16-26 0.888+ 2 

- 2 


16-26 0.728+ 4 

- 4 

16-26 0.805+ 3 

- 3 

16-26 0.849+ 3 

- 3 

18-26 0.901+ 2 

- 2 


16-26 0.810+ 4 

- 4 

16-25 0.861+ 3 

- 3 

16-26 0.886+ 4 

- 4 

16-25 0.924+ 3 

- 3 


16-26 0.877+ 4 

- 4 

18-26 0.906+ 5 

- 5 

18-26 0.920+ 6 

- 6 

18-26 0.946+ 5 

- 5 


17-26 0.945+ 7 

- 7 

17-26 0.961+ 6 

- 6 

17-26 0.966+ 6 

- 6 

17-26 0.979+ 6 

- 6 


16-26 0.986+ 5 

- 4 



TABLE VII: The table lists the masses of the lowest-lying, positive-parity spin-3/2 decuplet baryons at each value of the light 
(u/d) quark mass, together with the time range used in the fit. The quoted errors are obtained from the 68% confidence level 
using 10,000 bootstrap sub-samples. 









asqtad 

am , 

u/d 






Channel 


0.007 


0.010 


0.020 


0.030 


0.040 


0.050 


A 


16-22 0.948+ 15 

- 14 


17-22 0.998+ 16 

- 16 


16-26 1.007+ 9 

- 9 


16-26 1.058+ 10 

- 10 


16-26 1.072+ 12 

- n 


16-26 1.115+ 8 

- 8 


E* 


16-26 0.991+ 7 

- 7 


17-26 1.035+ 8 

- 8 


16-26 1.044+ 7 

- 7 


16-26 1.081+ 9 

- 9 


16-26 1.085+ 11 
- n 






16-26 1.047+ 4 

- 4 


17-26 1.069+ 6 

- 6 t 


16-26 1.082+ 6 

- 6 


16-26 1.104+ 8 

- 8 


16-26 1.099+ 11 

- 10 




0, 


16-26 1.101+ 3 

- 3 


18-26 1.113+ 4 

- 4 


16-26 1.119+ 5 

- 5 


18-26 1.118+ 10 

- 10 


16-26 1.112+ 10 

- 9 






IV. 


BARYON MASS CHIRAL EXTRAPOLATIONS 





In this section, we perform chiral extrapolations of the ground-state octet and decuplet baryon masses. We begin 
with a brief summary of our main results before presenting the details of our chiral extrapolation analysis in Sees. IV B 
and IV C. We use a generalization of heavy baryon chiral perturbation theory (HBxPT) [42, 65, 66] that includes the 
relevant lattice spacing effects for this mixed action lattice calculation [31]. We also perform a two-flavor extrapolation 
of the nucleon and delta masses, using continuum HBxPT. At 0(m 4 ), even the continuum two-flavor formulae have 
too many unknown low energy constants (LECs) to be determined from our lattice results alone. Therefore we must 
input various values to the formulae, including the nucleon, delta and nuclcon-dclta axial couplings, denoted in this 
work as gA, <?aa and <?A7V respectively. At this order, the two-flavor mixed action formula for the nucleon mass 
has an additional three unknown parameters, not including those from lattice spacing effects [38, 67], making the fit 
impossible with our results, even with the input of the LECs as in the continuum fit. Additionally, we only have 
numerical results at one lattice spacing, and therefore we can not control the continuum extrapolation. Despite these 
limitations, as explained below, we find at a pion mass of = 137 MeV (and itik = 497.6 MeV), 

M N = 954 ± 42 ± 20 MeV from NNLO SU{2) HBxPT 

= 960 ± 24 ± 8 MeV from NLO SU(6\S) mixed action HBxPT . (21) 

where the first uncertainty is statistical and the second uncertainty is an estimate of the systematic error in the 
extrapolation. While our extrapolated nucleon mass agrees with the physical value, using either the two-flavor 
continuum formula or the mixed action three-flavor formula, we must add a cautionary note. In the case of the 
three-flavor extrapolations, the values of the octet axial couplings, D and F, as well as the value of the octet-decuplet 
axial coupling, C, which we determine in the mass fits, arc significantly different from the known phenomcnological 
values fit with SU(3) HBxPT [68-73] and they are also inconsistent with the more recent lattice calculation of the 
octet axial charges [24]. The discrepancy with respect to the phenomenologically determined values of D, F and C 
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is independent of whether one uses the continuum SU (3) HBxPT formulae, or some generalization thereof, 4 raising 
questions about the convergence of the three flavor heavy baryon theory. 

A salient result of our calculations is the observation that, without any known theoretical explanation, the octet 
baryon masses are surprisingly linear in to^ over the range of pion masses studied in this work. For example, using 
the fit ansatz, 5 

M N = + cti m v , 

-a^ + a^^Tq, (22) 

and fitting to the lightest five of our nucleon mass results, we find 

M N = 938 ± 9 McV . (23) 

Furthermore, examining the nucleon mass calculations of several other lattice groups, we find that this is not unique 
to our results, and in fact this trend of the nucleon mass data is common to most other lattice calculations, employing 
a variety of different lattice actions, with different lattice spacings, some with three light dynamical quarks and some 
with two, and all with relatively light pion masses. 

Before presenting the details of our analysis, we first summarize the mixed action extension of HB%PT we use to 
perform the chiral extrapolations as well as discuss our general fitting procedure. Fitting is an issue with the three- 
flavor extrapolations because our analysis involves fits with up to 7 unknown parameters that are highly correlated. 

In heavy baryon x?T [42, 65], the octet baryon masses have a quark mass expansion given by 

M B = M + SM B 1] + SM { B 3/2) + SM B 2) + ... (24) 

where Mq is the quark mass independent contribution to the baryon masses in the chiral limit. The corrections, 
5M B n \ denote corrections to the mass of the baryon B, which scale as to™. The inclusion of the spin-3/2 states 
(resonances) in the theory somewhat complicates the expansion, since their inclusion introduces a new scale into the 
theory, the decuplet-octet (delta- nucleon) mass splitting in the chiral limit, which is a chiral singlet and leads to a 
renormalization of all parameters in the Lagrangian [67, 80, 81]. With these states included, one typically adopts a 
power counting [42, 65, 82] that includes the mass splitting 

A = Ma — M N ~ m n , (25) 

which is certainly relevant for our lattice calculations as well as the physical point. The dccuplet baryon masses then 
have a similar expansion given by 

M T = M + A + SM^ 1] + 5M^ /2) + 5M^ ] + ... (26) 

with A = A| m — o- At finite lattice spacing, these extrapolation formulae are modified in a known way. The two- 
flavor mixed action extrapolation formula for the nucleon has been determined to next-to-leading order (NLO) in 
Ref. [31]. To determine the formulae for the delta masses as well as the rest of the members of the octet and decuplet 
baryons, one can either extend the work of Ref. [31] or use the prescription in Ref. [38] to convert the known partially 
quenched formulae [67, 74, 80, 81, 83] to the relevant mixed action formulae. At NLO, one needs knowledge of 
the staggered taste-identity meson mass splitting [61] as well as the mixed meson mass splitting, which has been 
determined in Ref. [36]. Furthermore, supplemented by a treatment of the one loop flavor disconnected diagrams to 
the baryon masses [84], one can determine the mixed action extrapolation formula to next-to- next-to-leading order 
(NNLO), 0(m^) ~ 0(et 2 mj), in the mixed action effective theory. Here and throughout, we use a shorthand for the 
dimensionless expansion parameter, e 2 ~ ft2 A 2 ~ m 2 /A 2 , in which we drop the relevant A x 's. 



4 These include partially quenched [74] or mixed action [31] HBxPT, whether or not the decuplet degrees of freedom are explicitly 
included in the theory, and the so called "covariant" baryon x?T [75] as well as a different regularization scheme such as finite range 
regularization [76-79] . 

5 Plotting our resulting pion masses (squared) versus the input domain- wall quark masses, one can verify the expected scaling, ~ m q . 
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A. Fitting method for global chiral extrapolations 



When performing the chiral extrapolations, in addition to a standard \ 2 minimization, we perform, for lack of a 
better name, a bootstrap chiral extrapolation. As discussed above, for a given baryon mass, there are too many LECs 
in the NLO chiral extrapolation formulae to be determined from our lattice results alone. However, by considering a 
global x 2 minimization of all four non-degenerate octet baryon masses (or the four non-degenerate decuplet baryon 
masses), we can perform a fit to all seven (five) LECs that appear in the three flavor extrapolation formulae for 
the octet (decuplet) masses at NLO. In principle, we can even perform the NNLO fits that have 18 (14) LECs for 
the octet (decuplet) masses, although we find these fits generally do not converge and are not stable to changes in 
the fit ranges and the values of possibly fixed parameters. When performing the global \ 2 minimization, there are 
two correlations that are important to include, one of which makes these fits quite challenging. The first correlation 
is among the lattice results. For a given quark mass ensemble, m007, mOlO, etc., the various baryon masses are 
all correlated. Therefore one either should construct an error correlation matrix that accounts for this, or use the 
jackknife or bootstrap analysis techniques to handle these correlations. The second correlation is among the LECs 
being determined in the minimization, which is in fact the most challenging part of the analysis. In terms of the 
re-scaled dimensionless error correlation matrix, many of the off-diagonal elements are close to unity. Here we describe 
our bootstrap chiral extrapolation method that takes into account both of these correlations, and as far as we are 
aware, has not been implemented before. 

We begin by constructing bootstrap lists for all the baryon masses, with an equal number of bootstrap samples on 
each of the different quark mass ensembles, m007, mOlO, m020, m030, to040 and m050. For this work, we generated 
Ni, s = 10, 000 bootstrap samples at each quark mass. With these bootstrap lists, we then construct a bootstrap list 
of the global x 2 , 



where the sums run over B = {N, A, E, S} (and similarly for the decuplet), m q = {0.007,0.010,0.020,0.030, . . .}, up 
to the heaviest mass set used in a given minimization and the g(m,B ■ /irL m g]; m w [m q ], mK[m q ], . . . , A^) are the various 
baryon mass, m_e, chiral extrapolation functions that depend upon the masses and decay constants as well as the 
LECs, \i. In constructing the bootstrap lists of these global x 2 [Mj we nx the meson masses and decay constants, 
/j[m ? ], m,[mj and [m g ] to their central values, since their bootstrap fluctuations are an order of magnitude 
smaller than those of the baryon masses. We also do ignore the bootstrap fluctuations of the error of a given baryon 
mass, o- mB [m q ], taking their central values from tables II, III, VI and VII respectively. 6 To construct this list of 
global x 2 functions, as can be seen with Eq. (27), it is essential to have an equal number of bootstrap samples on the 
different quark mass ensembles, precluding the use of the jackknife method, at least with our sets of ensembles. 

The advantages of this method then follow naturally. Because the baryon masses on the different quark mass 
ensembles arc statistically independent, this amounts to adding independent noise, weighted by the statistical error 
fms [ m q]i to the central value of a given baryon mass on each of the different quark mass ensembles, mimicking a well 
known method of handling fits with highly correlated parameters. We note that for our calculations, the bootstrap 
samples appear Gaussian distributed about the mean. Furthermore, the correlations among the different baryon 
masses on a given quark mass ensemble are automatically taken into account by use of the bootstrap distributions. 
Therefore, by minimizing each of the Nb s entries in the global x 2 [M list; we generate a bootstrap list of the fit 
parameters, or the determined LECs, {A.;[frs]}. 7 Using the bootstrap error analysis, we can then make predictions for 
the resulting LECs, as well as the extrapolated baryon masses, which accounts both for the correlations among the 
ensembles at a given quark mass as well the correlations among the LECs from the minimization procedure, 



m B [bs, m q ] - g{m B : fn[mg],m„[m q ],mK[m g ], . . . , Xj) 

0~ms [^<?] 



2 



X 



2 n = EE 



(27) 



B m q 





(28) 



The inclusion of these fluctuations amounts to an error on the error, which is beyond our consideration here. 

We could also generate the bootstrap list of the error correlation matrix, but this would also amount to an error on the error, so we 
only retain the central values of this matrix. 
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and 

m B = m Bo ± a mB , m B() = g(m B : fv hys , m^ hys , m p ^ ys , X io ), 



\ A^TT E (™*o - ff(mn : ml hys , m% ys , A, [6s])) ' . 



To perform the minimization, we use both Mathematica and MINUIT. 

B. Two flavor chiral extrapolations 

In this section, we perform two-flavor chiral extrapolations of our nucleon and delta mass results. From the point 
of view of testing predictions from HB^PT, i.e. looking for non-analytic chiral behavior, one would ultimately like to 
determine values of qa-, 9aa and gAN directly from the nucleon and delta mass extrapolations. These LECs represent 
the leading order (LO) axial charges of the nucleon, delta and nucleon-delta transitions. The leading virtual pion 
cloud contributions to the nucleon and delta masses are proportional to these couplings, contributing at NLO. For 
example, the nucleon mass takes the following form at NLO 

M N = M - 2a M (jM) ml - tPt^^ ~ J/V^ Kr, A ' A > ( 30 ) 
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However, fits to our lattice results for the nucleon mass with this NLO formula return values of gA and gAN that are 
inconsistent with our knowledge of these LECs from either phenomenology or from lattice QCD. 8 In large part, this 
can be understood from the observation that our lattice data of the nucleon mass are well approximated by Eq. (22), 
Mjv = Q-o + c*i (see Table VIII and Figure 10 for more details). Therefore, in order for the SU(2) HBxPT 
nucleon mass expression to fit our lattice results, the different orders in the heavy baryon expansion of M/v, which 
is a polynomial scries in beginning at 0(?n 2 ) and supplemented by chiral logarithms, must conspire to form this 
straight line. Viewed in this light, it is not surprising that an NLO analysis fails to return values of gA and gAN that 
are consistent with their known values. Therefore, in accord with the expectations from Refs. [85, 86], we must use 
the extrapolation formula to at least NNLO and ideally next-to-next-to-next-to-leading order (NNNLO) to test if the 
values have stabilized. At this order, unfortunately, there are too many LECs in the formula to be determined from 
a fit to the lattice data points alone, a problem that is only exacerbated with the use of the mixed action formula. 
Hence, we must resort to fixing the values of some of these parameters using results either from phenomenology or 
other lattice calculations, ideally determined with the same lattice action. 

1. The nucleon mass 

To study the chiral extrapolation of the nucleon mass, we perform fits to our lattice data with the following formulae; 

• LO: SU{2) HBxPT, 

M N = Mq - 2a M ml , (32) 

• NLO: SU(2) HB^PT, Eq. (30) [42, 65, 66]. 



The use of the mixed action expression for the nucleon mass at this order [31], supplemented by the known valence-sea meson mass 
splitting [36] does not qualitatively change this conclusion. See Table VIII for details. 
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NLO: SU(A\2) mixed action HB X PT [31] with 
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where gi is an additional axial coupling appearing in the mixed action/partially quenched Lagrangian [83], 
Apg = a 2 Ai is the taste-singlet staggered meson splitting, and m 2 u = m 2 + a 2 AMix, with the mixed valence- 
sea meson mass splitting given in Ref . [36] . 

• NNLO: ST/ (2) HB^PT with explicit delta degrees of freedom (explicit deltas) [67] (and a slightly modified 
renormalization prescription from Ref. [67] ) 9 
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• NNLO: SU{2) HBxPT without explicit deltas, Eq. (34) with g AN = 0. 

• NNLO: SU(2) covariant baryon x?T formula without explicit deltas, expanded to 0(m%) [87] 
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• Mpf = Q!q + oti m^: an empirical form motivated by the observed nucleon mass results, not motivated by any 
understanding of low energy QCD we currently have. 

Performing the LO through NNLO fits in principle allows us to study not only the chiral convergence of the nucleon 
mass, but also to monitor the resulting values of the LECs as higher order terms are added to the expansion. If 
things are working as desired, we would find not only that the expression for the nucleon mass is converging order 
by order, but that the values of the LECs determined in the analysis would shift by only small amounts as we add 
higher order terms. Unfortunately, since we must fix certain LECs in the NNLO fits, we are not able to honestly make 
this comparison. In the NLO and NNLO expressions, the LECs f^, gA and gAN can either take their LO values, 
or one can replace them with their lattice-physical values, 10 the difference appearing at NNNLO. To maintain as 
much consistency between the various fits as possible, including analyses performed by other groups, we always take 
/tt = 121.9(8.8) MeV, consistent with the two loop determination of f n in the chiral limit [88]. We consider values of 
gA that are consistent with the physical value, the phenomenological value in the chiral limit, and the lattice value of 



The coefficients b^(ii) and &A:f(/-0 are not the renormalized coefficients as defined in Ref. [67]. There are additional operators which 
contribute to the nucleon mass at NNLO, notably the operators which are responsible for the LO contribution to the delta mass. 
However, the contribution to the nucleon mass from these other operators is parametrically (in m^) the same as those proportional to 
bA and 6m- Since we are not doing a combined fit of the nucleon and delta mass, we have absorbed these effects with a re-definition of 
these coefficients. 

We denote lattice-physical quantities as those that are determined directly from correlation functions, and that have not been extrapo- 
lated to the continuum, infinite volume or physical/chiral point. 
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TABLE VIII: Various two-flavor nucleon mass chiral extrapolation results. In the left column, we provide the particular fit 
function as well as the values of LECs fixed in the minimization. In the far right column, we provide the resulting nucleon 
mass at — 137 MeV. In the fit parameters as well as the predicted nucleon masses, the first uncertainty is statistical. In the 
NNLO fits, the second uncertainty is systematic determined by varying the fixed LECs over their given ranges. The last fit 
function is motivated purely by the observed lattice results for the nucleon mass, and not by any understanding of low-energy 
QCD we currently have. We have set the renormalization scale to /i = 1 GeV. 



FIT: LO 


range 


Mo [GeV] 


car [GeV -1 ] 






X 


d.o.f. 


Mjv [MeV] 


Mjv = Mo - 2a M ml 


007-020 


1.00(1) 


-0.57(3) 






1.4 


1 


1028 ± 9 




007-030 


1.02(1) 


-0.53(2) 






4.6 


2 


1037 ± 8 




007-040 


1.02(1) 


-0.51(1) 






6.8 


3 


1043 ± 7 




007-050 


1.04(1) 


-0.47(1) 
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4 


1056 ± 6 


NLO SU(2), Eq. (30) 


range 


Mo [GeV] 


a M [GeV L ] 


9A 


g a at 


2 

X 


d.o.f. 


M at [MeV] 


fo = 121.9(8.8) MeV 


007-040 


0.98(2) 


-0.80(12) 


0.43(9) 


0.00(1.86) 


1.39 


1 


1013 ±15 




007-050 


0.98(1) 


-0.84(8) 


0.47(6) 


0.00(2.48) 


1.60 


2 


1009 ± 12 


NLO 5(7(412), Eq. (33) 


range 


Mo [GeV] 


a M [GeV -1 ] 




g^N 


x 2 


d.o.f. 


Mjv [MeV] 


fo = 121.9(8.8) MeV 


007-050 


1.01(4) 


-0.95(13) 


(0.6(2) , -0.5(1.4)) 


0.03(3.60) 


1.66 


1 


1046 ± 38 


NNLO, Eq. (34) 


range 


Mo [GeV] 


ctM [GeV -1 ] 


b M [GeV- 3 ] 


b A 


x 2 


d.o.f. 


Afjv[MeV] 


fo = 121.9(8.8) MeV 


007-040 0.87(6)(3) 


-3.1(7)(8) 


62(11)(30) 


-29(8) (16) 


0.06 


1 


941 ± 42 ± 17 


g A = 1.2(1), g AN = 1.5(3) 


007-050 0.90(4) (5) 


-2.7(4)(9) 


55(7) (24) 


-24(5)(17) 


0.75 


2 


966 ± 43 ± 20 



NNLO, Eq. (34) 


range M [GeV] 


olm [GeV -1 ] 


6m [GeV -3 ] 


bA X 2 


d.o.f. 


Mjv [MeV] 


fo = 121.9(8.8) MeV 


007-040 0.91(5)(0) 


-1.8(5)(1) 


4.6(0.4)(1.0) 


-7.6(4.2) (0.9) 0.00 


1 


964 ±41 ±20 


gA = 1-2(1), gAN — 


007-050 0.96(4) (5) 


-1.4(3)(6) 


4.8(0.3)(1.1) 


-3.6(2.3) (4.1) 1.36 


2 


996 ± 30 ± 30 


NNLO, Eq. (36) 


range Mo [GeV] 


Cl [GeV -1 ] 


ei[GeV -3 ] 


x 2 


d.o.f. 


Mjv [MeV] 


fo = 121.9(8.8) MeV 


007-030 0.90(2) (1) 


-0.97(4) (8) 


2.8(5)(9) 


0.02 


1 


958 ± 15 ± 9 


g A = 1.2(1) 


007-040 0.90(1) (2) 


-0.97(2) (8) 


2.7(2)(8) 


0.07 


2 


956 ± 12 ± 11 


[ C2 ,c 3 ] = [3.2(4), -3.4(4)] GeV" 


1 007-050 0.88(1)(2) 


-1.01(1)(10) 


2-2(l)(8) 


0.5 


3 


940 ± 9 ± 14 



Eq. (22), Mjv = oto + af range 


a^GeV] 


«i 


x 2 


d.o.f. 


Mjv [MeV] 


007-020 


0.83(2) 


0.93(5) 


0.00 


1 


953 ± 13 


007-030 


0.82(2) 


0.94(4) 


0.18 


2 


950 ±11 


007-040 


0.80(1) 


0.99(3) 


4.39 


3 


938 ±9 


007-050 


0.80(1) 


1.01(2) 


5.40 


4 


933 ±8 



the nucleoli axial charge calculated with this mixed action approach [5], gA = 1.2(1). see also Ref. [89]. To fix gAN-, 
we use the known width of the delta, combined with the LO expression for the width from HBxPT [65, 67] 

T A = -2Im[M A ] = |% ( A 2 - mlf' 2 , (37) 

from which we obtain <?ajv = 1.5(3). 11 The error we assign also encompasses the NLO determination of g^N [69]. We 
also set A = 271 MeV as well as insert the lattice value for A = MA[m g ] — M/v[m,j] into our extrapolation analysis 
for each quark mass. When using the NNLO formula, Eq. (36), we follow Ref. [92] and set c 2 = 3.2(4) GeV -1 and 
C3 = —3.4(4) GeV -1 (we have slightly inflated the errors from those in Ref. [92]). 



11 There are various sign conventions for the axial couplings used in the literature [90, 91]. The nucleon and delta masses are proportional 
to the square of these couplings, and we take them to be positive. 
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FIG. 10: Chiral extrapolations of the nucleon mass, plotted vs the approximate expansion parameter m^/2V27l'/o, with 
fo = 121.9 MeV. For comparison purposes, in all figures, we display the results of fits to the m007-m040 mass points denoted 
by the small black circles with error bars. The black square is the m050 mass point not included in any of these fits. The 
filled (red) circle is the physical nucleon mass, taken to be 939.6 MeV, at a pion mass m w = 137 MeV, which is never included 
in the minimization. The gray bands represent the 68% confidence interval, and are only determined from the statistical error 
bar in the lattice results. In Fig. 10(a), we plot both the LO and NLO SU(2) HBxPT results. The light shaded band is from 
LO and the darker shaded band is NLO. In Fig. 10(&), we plot the results of the NNLO 57/(2) fit including explicit deltas. 
In Fig. 10(c), we plot the NNLO 577(2) covariant fit without deltas. In comparing Figs. 10(6) and (c) one needs to note the 
size of the error band is dictated by the number of free parameters and not by the use of infrared-regularization (covariant 
expression). In Fig. 10(d), we plot the straight line fit, Eq. (22). All of these fits, except the LO fit of (a), are statistically 
consistent with our lattice results, as can be seen in Table VIII. 



We collect the results of these various extrapolations in Table VIII. The first uncertainty in the fit parameters 
and the predicted nucleon mass is statistical. In the NNLO fits, the second error is a systematic error obtained by 
independently varying the fixed LECs over their given ranges. In Figure 10, we display some of the resulting fits along 
with their statistical 68% confidence bands. For comparison purposes, we have chosen to display fits that include the 
m007-m040 mass points. 12 



If we instead perform these fits in ri units, as commonly employed by the MILC Collaboration, the qualitative features do not change. 
For example, converting our nucleon and pion mass results to rj units using Table 1 of Ref. [2] , using the m007-m040 ensembles, we 
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*a=1.2(1), gAN=L5(3) ^=1.2(1), gAN=l-5(3) 




m n (GeV) (GeV) 

(a) (6) 

FIG. 11: Plot of the contributions to the nucleon mass order-by-order, SM^ resulting from the NNLO heavy baryon 
fit including the delta. In (a) we plot each order separately, plotting the magnitude of the (negative) NLO results denoted by 
a dashed line. In (b) we plot the combined LO+NLO and LO+NLO+NNLO results. Plotting any of the NNLO fits of the 
nucleon mass in this work results in similar values for the 5M^ . The arrows denote the values of the pion mass used in this 
work. From (a), one notes that the convergence of the resulting NNLO fit may already be breaking down at the second lightest 
mass point, although the summed contributions in (6) display better convergence. 



From the analysis presented in Table VIII, we draw the following conclusions. The LO nucleon mass formula, 
Eq. (32), does not provide a good description of our lattice data, except for the lightest three mass points. At NLO, 
we see that both the continuum and the mixed action formulae provide reasonable descriptions of our lattice data, 
such that with the available degrees of freedom, we can not strongly conclude one fit is better than the other. However, 
they return values of qa and <?ajv that are inconsistent with both phenomcnological and lattice determinations. At 
NNLO, with fixed axial couplings, we see that the explicit inclusion (gAN = 1-5(3)) or exclusion (gAN = 0) of the delta 
does not have a significant impact upon the quality of the fit, or the chirally extrapolated nucleon mass. However, as 
expected, the values of the LECs au and more notably 6m and bA are very sensitive to this change. The NNLO fit 
based on covariant baryon %PT with infrared regularization, without explicit deltas, in fact has only three unknown 
fit parameters, as compared to the other NNLO fits, and this is the reason the statistical error bars are much smaller. 
The reason we chose to fix more parameters in the covariant NNLO fit is to have a more direct comparison with the 
chiral extrapolation analysis of other groups, who also chose to fix the values of c-i and c 3 from phenomcnological 
determinations [92-94]. It is worth noting that the values of the fit parameters using the NNLO covariant fit agree 
within errors with those determined in Refs. [87, 92, 93] which performed the same fit to two-flavor lattice calculations. 

When we let either of these LECs float as a free parameter in the extrapolation, then the fit is comparable to 
the NNLO HB^PT analysis we have performed. In Figure 11(a), we plot the resulting contributions to the nucleon 
mass, order by order for the NNLO fit with explicit deltas and in Fig. 11(6) we plot the combined LO+NLO and 
LO+NLO+NNLO mass contributions. A similar plot with any of the HB^PT NNLO fits returns a similar expansion, 
see also Ref. [95]. For visual aid, the arrows denote the values of the pion masses used in our lattice calculation. As 
can be seen in Table VIII and Figs. 10 and 11, even though the resulting fits to the nucleon mass are in good statistical 
agreement with our lattice data, already at the lightest value of the pion mass, Fig. 11(a) suggests the ordcr-by-ordcr 
convergence is becoming questionable. In contrast, Fig. 11(6) displays much better convergence, which is understood 
from the oscillatory nature of the expansion. 

The last fit we perform is a linear fit in ~ \J m q-, Eq. (22). This ansatz is not motivated by any theoretical 
understanding of low energy QCD, but rather by the empirical observation that the results are nearly linear in m w . 
It is then natural to ask how good this fit is statistically. It does not describe our lattice data as well as the NNLO 
analysis, however purely from the x 2 j this fit can not be ruled out. Here we would like to stress that we are not 
advocating Eq. (22) as a means of performing the chiral extrapolation. This fit ansatz is conceptually incorrect 
near the chiral limit where chiral symmetry dictates the nucleon mass scales as Mjv = Mo + a'm q . Rather we are 
highlighting this fit to bring attention to this unexpected phenomenon, and to ask the question, for which values of 



find af = 1.03(3) as compared to the value in Table VIII. The NNLO SU{2) HBxPT analysis with g A = 1.2(1) and gAN = 1-5(3) also 
does not change qualitatively, returning for example Mq = 0.93(6) GeV and om = —2.5(3) GeV -1 . 
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FIG. 12: Nucleon mass results from other lattice groups, the MILC [96], QCDSF/UKQCD [93], RBC/UKQCD [97] and 
ETM [94] Collaborations. We apply our linear fit in to the lattice results from all groups and plot the corresponding best 
fit and 68% confidence band. In all cases, the open diamonds correspond to the points included in the analysis. In Fig. 12(a), 
in addition to the coarse MILC (a ~ 0.124 fm) nucleon mass results, we also display the same analysis from our lattice results 
for comparison, as well as the MILC superfine (a ~ 0.06 fm) results, which were not used in the analysis. In Fig. 12(d), we also 
display a fit to the ETM lattice data using Eq. (36) , which provides a better \ 2 description of their nucleon mass results. 



TOjr will this fit ansatz fail to describe the nucleon mass? 

It is important to point out again that our lattice calculations have been performed with only one lattice spacing and 
only one lattice volume. The volume effects arc expected to be about 1% or less, for our lightest mass as determined 
at NLO in heavy baryon xPT [98], and therefore do not impact this analysis significantly. However, we have no 
control over the lattice spacing effects, and since we can not perform the NNLO mixed action analysis to compare 
to the continuum NNLO analysis, we can not provide more than an order of magnitude estimate of their size. For 
these reasons, it is important for us to compare our nucleon mass results with other groups, which use a variety of 
different lattice actions, at different lattice spacings, with slightly different lattice volumes and with various scale 
setting procedures. 

Recently, four other groups have reported their computational results for the nucleon mass, and in most cases 
the resulting chiral extrapolations: the MILC [96], QCDSF/UKQCD [93], RBC/UKQCD [97] and most recently the 
ETM [94] collaborations, all of which have employed relatively light pion masses in their lattice calculations. The 
MILC collaboration uses Nt = 2 + 1 flavor staggered valence fermions on rooted, staggered sea fermions, for which the 
nucleon mass is expected to have complicated lattice spacing dependence [99]. The MILC Collaboration has results 
for three lattice spacings, and here we compare to the coarse (a ~ 0.124 fm) and superfine (a ~ 0.06 fm) MILC results 
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from Rcf. [96]. The QCDSF/UKQCD group results are from Nf = 2. O(a) improved Wilson fermions, with several 
different lattice spacings, ranging from a ~ 0.085 — 0.067 f m, for which the nucleon mass is expected to have very simple 
lattice spacing corrections [100]. The RBC/UKQCD group employed Nf = 2 + 1 domain- wall on domain- wall fermions 
with a ~ 0.114 fm, and are expected to have very simple and mild lattice spacing corrections. The ETM collaboration 
has used Nf = 2, O(a) improved twisted mass fermions, and reported results for three lattice spacings. The nucleon 
mass with twisted mass fermions at maximal twist is expected to have very simple lattice spacing dependence [101], 
but more importantly, they do not find any significant lattice spacing corrections in their results [94]. In Fig. 12 we 
plot the resulting nucleon mass calculations from these groups along with our own straight-line in analysis of 
their data. The only data for which the SU (2) nucleon mass extrapolation formula performs significantly better than 
Eq. (22) are those of the ETM collaboration. Even in this case, the straight line is still a reasonable approximation 
as measured by the x 2 . 

Within statistical error bars, our resulting nucleon mass calculations at a ~ 0.124 fm are consistent with those 
of the MILC superfine (a ~ 0.06 fm), RBC/UKQCD (a - 0.114 fm) and QCDSF/UKQCD (a ~ 0.085 - 0.067 fm) 
nucleon mass calculations, while the ETM nucleon mass results are systematically higher. Our results are significantly 
different from the MILC calculations on the coarse lattices on which we performed our calculations. The systematic 
difference with these coarse MILC results is highly suggestive of an expected overall additive 0(a 2 ) difference in the 
nucleon mass from the two different actions. Given the common method of scale setting between our calculation and 
that of MILC, and given the agreement of our results with the superfine MILC results, 13 the DWF valence on the 
coarse MILC sea action seems to have significantly smaller discretization effects than the Asqtad action. 

As is abundantly clear from this analysis, this linear behavior in is not unique to our mixed action calculation 
with domain-wall valence fermions on the MILC staggered sea configurations. Is QCD conspiring to produce this 
remarkably straight line? Or is this flattening of the chiral curvature of Mm a lattice artifact? Unfortunately, there 
are still too many variables to answer this question. Half of the groups have used smaller lattice spacings, but only 
two flavors of light quarks. The other half have used three dynamical light quarks, but employ larger lattice spacings 
(except for the superfine MILC results which have the smallest lattice spacing), albeit with very different lattice 
actions and expected lattice spacing dependence. It is clear from Figs. 10, 11 and 12 that lighter pion masses are 
desirable, to determine if the lighter point would favor more curvature in the extrapolation function, although the 
lightest MILC superfine result has m n ~ 220 MeV. Also, it will be important to use a second lattice spacing with our 
mixed action as well as to perform complete domain wall calculations with domain wall sea quarks to explore the 
continuum extrapolation. It is worth mentioning that each group has chosen an independent means to set the scale 
of their lattice calculations, so this phenomenon is not an accidental feature of a particular scale setting method. 
To resolve this issue, a careful study with multiple lattice spacings, volumes and light pion masses will be needed. 
Another feature that is now quite clear is that the nucleon mass is a poor observable to use to look for signatures 
of expected chiral non-analytic behavior, predicted from HBxPT, at least at the pion masses used in this work and 
other recent lattice calculations. 14 



2. The (pion) -nucleon sigma term, <tjv 



We comment briefly on the pion-nucleon sigma-term, defined as the nucleon scalar form factor at zero momentum 
transfer. From the Feynman-Hellmann theorem, we can relate this to the slope of the nucleon mass with respect to 
the quark mass, 



Up to NNLO differences, this is equal to 



d 

cr N = m q -^—M N . (38) 



5n = ^J-M n . (39) 

2 OTOjr 



While the explicit inclusion of the delta states has little impact upon the value of the chirally extrapolated nucleon 
mass, it has a significant effect upon the sigma-term. For example, using the resulting NNLO HBxPT fits, Eq. (34), 



We thank Doug Toussaint and other members of the MILC Collaboration for providing us with their preliminary superfine nucleon mass 
results, some of which are not yet published. 

For an updated comparison and analysis of this phenomenon, see Ref. [102] . 
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TABLE IX: Various two-flavor delta mass chiral extrapolation results. We provide the minimization results for the NNLO fit 
to the delta mass using Eq. (42), fitting directly to our lattice results as well as to the results with the predicted finite volume 
corrections subtracted, as discussed in the Sec. IV B 3. 



NNLO, Eq. (42) 


range A/ A ,o[GeV] 


7A/[GeV _1 ] 


MGeV" 3 ] 


t A 


x 2 


d.o.f. 


A/ A [MeV] 


fo = 121.9(8.8) MeV 


007-040 1.59(13)(1) 


-1.1(2.5)(0.7) 4.5(1.8)(4.5) 


-27(33) (6) 


0.34 


1 


1570 ± 92 ± 58 


9 an = 1-5(3), 3aa =2.2(6) 


007-050 1.49(9) (2) 


1.1(1.6)(0.8) 


6.4(9)(4.7) 


3(20)(2) 


1.66 


2 


1500 ± 92 ± 78 


NNLO, Eq. (42) with predicted 5$ v 


range A/ A ,o[GeV] 


7A/[GeV _1 ] 


t M [GeV- 3 ] 


t A 


x 2 


d.o.f. 


A/a [MeV] 


fo = 121.9(8.8) MeV 


007-040 1.52(15)(1) 


-1.0(2.9)(0.7) 4.4(2.1)(4.9) 


-28(39) (6) 


0.01 


1 


1509 ± 92 ± 6 


9 an = 1.5(3), g a a = 2.2(6) 


007-050 1.38(9)(10) 


2.0(1.8)(0.9) 


6.7(9)(4.8) 


13(22)(10) 


1.66 


2 


1412 ±92 ±9 




007 010 


020 


030 










S£ v = (Af A [2.5 fm] - M A [oo])[MeV] 


42 ± 18 39 ± 16 


4.9 ±2.2 


1.8 ±0.8 











wc find 



a N = 84 ±17 ±20 MeV, 
a N = 42 ± 14 ± 9 MeV, 



with gAN = 1-5(3) , 
with gAN = , 



(40) 



taking an average of the m007-m040 and m007-m050 fits. While these fits are consistent at the 68% confidence 
level, the explicit inclusion of the delta states increases the central value of the sigma-term by a factor of two. With 
gAN — 1-5(3), the systematic error is dominated by this uncertainty in gAN- Therefore, to really determine the 
nucleon sigma term from lattice QCD, one will need lattice data from which one can reliably determine gA and gAN 
from the chiral extrapolation of the nucleon mass, or alternatively, lattice results in close proximity to the physical 
point on both sides. Just for comparison purposes, we find, using the straight-line fit, a value of aN = 67 ± 2 MeV. 



3. The delta mass 



The chiral expansion of the delta mass takes a similar form to the nucleon mass, 
extrapolation formula for the delta mass is given by [67] 



At NLO, the SU{2) HB^PT 



M A = M + A + j M ml - 
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(41) 



For < A, !F(m w , —A, fx) develops a branch cut, producing the imaginary part of the delta self-energy responsible 
for the strong decay to a p-wave nucleon-pion state. In the finite box we work with, where the particle momenta 
are quantized in units of p = 2tt/L ~ 500 MeV, this decay is kinematically forbidden. In the center of mass frame, 
the nucleon and pion must each have about p ~ 266 MeV, and one might expect that as in the case of a delta that 
would be stable in infinite volume, the finite volume corrections to the delta mass should be exponential, exp(— m^L). 
However, as noted in Ref. [103], since the delta undergoes a strong decay and its width is therefore large, the plateau 
regions of the energy levels vs. L characterizing narrower resonances [58] are "washed out" for the delta. According to 
Ref. [103], which uses an NLO covariant baryon xPT calculation, the energy levels receive power law in L corrections, 
rather than exponential. In this work, we perform two different extrapolations. First, we perform the fit with the 
continuum, infinite volume extrapolation formula. Second, wc subtract from our lattice results the finite volume 
corrections calculated according to Ref. [103], and then perform the chiral extrapolations. Ultimately, a careful study 
with multiple volumes is required to test these predicted corrections and perform reliable extrapolations of the delta 
mass. 

We note that recently some groups have opted to perform the chiral extrapolation of the delta mass by setting 
the strong coupling of the delta to the nucleon-pion state to zero, by hand, accomplished by setting gAN = in 
Eq. (41) [93, 94]. We chose not to follow this approach. To perform the chiral extrapolation of the delta mass, we 
use the continuum SU(2) HBxPT extrapolation formula to NNLO [67], with a similarly modified renormalization 
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TABLE X: Resulting delta mass determinations on the 20 3 x 64 and 28 3 x 64 coarse MILC ensembles, using both the smeared- 
smeared (SS) as well as smeared-point (SP) correlation functions. The results on the unchopped 20 3 x 64 ensembles are in 
agreement with our high statistics results on the chopped 20 3 x 32 ensembles. Note the source on the 20 3 x 32 ensembles was 
placed at t = 10. 



SS 



Volume 
20 3 x 32 
20 3 x 64 
28 3 x 64 



Range 
16-22 
1-10 
1-10 



0.974(9) 
0.95(3) 
1.01(4) 



Range 



10 
10 



SP 



0.97(3) 
1.03(3) 



prescription as discussed for the nucleoli mass. Taking the real part of the mass formula, we obtain 
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(42) 



(43) 



(44) 



As with the nucleon mass, there are too many LECs to be determined from the lattice results alone. We therefore fix 
9 an = 1-5(3) as above and c/aa = 2.2(6) [69-71]. When performing our delta mass extrapolations, we take the results 
of the fits that include oscillatory terms. Comparing Tables V and VII, one can see that the resulting delta mass on 
the mOlO ensemble is smaller using Eq. (10) than with the single exponential fit (this discrepancy is negligible for all 
other ensembles). Examining Fig. 9, it is clear that the oscillatory region provides constraints with small error bars, 
so that a fit including the oscillatory terms discussed in Sec. Ill A provides a better description of the correlation 
function, which is confirmed by the x 2 minimization. We collect the results of our analysis in Table IX, as well as the 
predicted finite volume corrections, defined as 



8f v ee M A [2.5 fm] - M A [oo] 



(45) 



These values have been obtained using the central values of m^, Mjv and Ma, listed in Tables II, VI and V respectively 
and by varying gAA, 9 an and / w within the ranges specified in the present paper. In Figure 13 we present the resulting 
chiral extrapolations. In the figure, the open circles with error bars arc our lattice results. The open squares that sit 
slightly below the circles are the lattice results after subtracting the predicted finite volume corrections, determined 
as described in Ref. [103]. The filled circle with an error bar is our result for the delta mass using the mOlO ensembles 
with L ~ 3.5 f m lattices, for which we display the correlation functions and oscillating fits in Fig. 14, and the resulting 
masses in Table X. Note, this lies in the opposite direction as that predicted from Ref. [103]. However, the minimum 
momenta on these lattices are p ~ 350 MeV. The pion-nucleon p-wave decay channel is therefore more accessible and 
this state may have much larger volume corrections. Clearly, a multiple volume study will be required to study the 
delta mass and its chiral extrapolation. 15 



A fit using the NNLO covariant baryon xPT delta mass extrapolation formula, as advocated in Ref. [104], might be in more agreement 
with the pole position of the delta mass. However, given our numerical results of the delta mass in the two different volumes, there are 
clearly more important systematic effects to be understood than which version of baryon EFT provides a better agreement with the 
physical pole mass. We therefore do not make use here of the formula in Ref. [104], 
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FIG. 13: Chiral extrapolations ol the delta mass. In (a) we display the extrapolation of our lattice data using the continuum 
infinite volume formula. In (b), we display the extrapolation of our lattice data with the predicted finite volume corrections 
subtracted from our lattice results. For comparison, both fits are done using all six mass points. The open circles with error 
bars are our calculational results while the lower open squares are the predicted infinite volume extrapolations. The filled 
circle with error bar represents our preliminary determination of the delta mass using the mOlO ensemble on the L ~ 3.5 fm 
lattices. Note that it is located on the opposite side compared to the predicted infinite volume extrapolation of Ref. [103] . The 
filled (red) circle denotes the pole mass of the delta plotted at m, = 137 MeV. A multiple volume study is required for the 
extrapolation of the delta mass, which is beyond the scope of this work. 



Lattice size: 20 

Smearcd-smcared 



Lattice size: 28~ 

Smcared-smeared 





(a) 



(b) 



FIG. 14: Oscillating fits to the smeared-smeared delta correlation function, Eq. (10) on the mOlO 20 3 x 64, (a) and the 28 3 x 64, 
(b) lattices. The delta mass on the larger volume is heavier than on the lighter volume, contrary to expectations. 
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C. Three flavor chiral extrapolations 



We now turn to the chiral extrapolations of the octet and decuplet baryon masses, making use of the mixed 
action generalization of SU(3) heavy baryon xPT, as discussed above. 16 Before delving into the details of our chiral 
extrapolations, we first discuss our strategy. In addition to fitting the baryon masses, we also fit the baryon mass 
splittings, both exclusively and in combination with the masses. There are several advantages in studying the baryon 
mass splittings vs. the masses, a technique first successfully exploited in Ref. [107] in a quenched baryon spectroscopy 
calculation. First, the leading lattice spacing corrections to the masses exactly cancel in the splittings, pushing these 
lattice artifacts to NNLO in the mixed action EFT [31, 38]. 17 Second, the mass splittings must identically vanish 
both in the three-flavor chiral limit as well as the SU(3) degenerate mass point. This is reflected in the extrapolation 
formulae and leads to smaller theoretical error in the extrapolation to the physical point, as has been successfully 
employed for other quantities with similar limits calculated with chirally symmetric valence fermions [12, 16, 17, 19, 20]. 
Finally, fitting the mass splittings as well as the masses allows us to address in more detail the convergence of SU (3) 
heavy baryon x?T, a topic that we address after detailing the mixed action Lagrangian. 

The partially-qucnchcd/mixed-action Lagrangian in the notation of Refs. [38, 74], is given by 
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(46) 



The last term in each of the first two lines provides the leading lattice spacing correction to the baryon masses. As 
can be seen, these terms treat all the octet and decuplet baryon masses the same, and therefore, as mentioned above, 
this leading lattice spacing dependence drops out of the mass splittings. For our purposes, we only need the leading 
term of the mass spurion field, 



M. 



(47) 



where the dots represent terms with higher powers of the meson fields, and m q is the quark mass matrix. Here, we 
have suppressed all flavor indices, and the braces ( ) represent the flavor traces as defined in Ref. [74], and str() stands 
for a super-trace over the graded flavor algebra. The relation between the LECs of this Lagrangian and the standard 
SU(3) Lagrangian [42, 65] can be determined by matching this theory onto the valence sector [74], for which one finds 
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In Ref. [105], the chiral extrapolations of the octet baryon masses calculated by MILC on the coarse Asqtad ensembles was explored 
using the NNLO covariant formula determined in Ref. [106] . 

The virtual meson cloud contribution to the baryon masses, which enters the mass expressions at NLO, receives some contributions from 
mixed- mesons that have discretization errors in them. However, the additive mixed meson mass splitting has been determined [36], so 
this lattice artifact introduces no new unknown LECs, as opposed to those that enter at NNLO. 
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From this Lagrangian, for example, the nucleoli mass is given to NLO by [74] 
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where T{m^ 1 A, /i) is defined in Eq. (31) and T(m$, 0, /i) = irm^. We use modified LECs 
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that have mass-dimension —1, and also we subsequently drop the primes. The parameter B$ is related to the chiral 
condensate that appears in the Gell-Mann-Oakes-Renner relation [108] and the meson chiral Lagrangian. We further 
have 
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where a 2 Ai and a 2 AMix can be found in Ref. [61] and [36] respectively. Similar expressions for all the octet and 
decuplet baryon masses are determined using Refs. [31, 67, 74, 80, 81, 83, 84]. We use these formulae to perform 
the chiral extrapolation analysis of our baryon mass results, that we turn to after a brief digression on the expected 
convergence of SU(S) HBxPT. 



1. On the convergence of SU(3) heavy baryon xPT 
The expansion parameter of xPT [39-41] in the meson sector is given in terms of generic meson masses in the set 

{?7l0 : 777^,777^,777,,}, aS 
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~ -A ■ (52) 

x 

The inclusion of heavy matter fields in the effective theory, such as baryons or heavy mesons, leads to the expansion 
parameter 

e H ~^, (53) 

for which the convergence of the theory becomes worse and/or questionable. In fact, using general knowledge of 
asymptotic scries and taking the ratio of physical parameters, m n /A x ~ 1/2, one may expect that SU(3) heavy 
baryon x?T would only have a chance of converging for 2-3 orders before the asymptotic nature of the theory became 
apparent. This problem is exacerbated with current lattice calculations, in which the meson masses are heavier than 
those in nature, casting doubt on the ability to reliably extrapolate observable quantities computed on the lattice, 
using three-flavor heavy baryon xPT, or any of its generalizations. 18 



Those generalizations include explicit inclusion of the decuplet degrees of freedom, partially quenched or mixed action generalizations 
for the lattice, covariant baryon xPT with infrared regularization [75] , or some other regulator like finite range regularization [76-79] . 
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TABLE XI: NLO, SMg^ 2 \ contributions to octet baryon masses and mass splittings in heavy baryon xPT both with decuplet 
degrees of freedom and without We use the charge-neutral baryons in these relations as the lattice computations without 
electromagnetism are most similar to these. To determine the predicted NLO mass contributions and their corresponding 



errors, we use D = 


0.715(50), F = 


= 0.453(50) and C = 


1.2(2). 








Quantity 


Experimental 


HBxPT 


4 (c = o) 


Quantity 


Experimental 


HByPT A (C = 0) 






6m£/ 2) - SM% /2) 


SMg/ 2 - 1 - SM ( B 3/2) 






*Mj, 8/3) 5Mi 3/2) 


M A - M N (MeV) 


176 


-285(65) 


-326(65) 


Mjv (MeV) 


940 


-195(38) -278(38) 


Ms - M N (MeV) 


253 


-152(60) 


-287(46) 


Ma (MeV) 


1116 


-480(84) -604(84) 


Ms - M N (MeV) 


375 


-516(120) 


-637(120) 


Me (MeV) 


1193 


-347(79) -565(66) 


M E - M A (MeV) 


77 


133(76) 


39(76) 


M H (MeV) 


1315 


-711(124) -915(124) 


Ms - Ma (MeV) 


199 


-232(39) 


-311(39) 








Mb - Me (MeV) 


122 


-367(81) 


-351(81) 








Aqmo (MeV) 


10 


9(4) 


3(4) 









Recently, the first calculation of the hyperon axial couplings from the lattice was performed, in which it was found 
that 5J7(3) heavy baryon xPT at NLO failed to describe the lattice results [24]. Additionally, the RBC/UKQCD 
Collaborations have advocated the use of two-flavor chiral extrapolations for their recent calculations of meson quan- 
tities including the decay constants as well as the kaon bag parameter [109]. In this work, we find that both SU(3) 
heavy baryon xPT as well as its mixed action generalization describe the lattice results of the baryon masses and 
mass splittings, with reasonable \ 2 per d.o.f. However, these fits (presented in Sec. IV C 2) return values of the axial 
coupling LECs, D, F and C that are far from the recent lattice calculations [9, 24] as well as the phenomcnologically 
determined values [69-73], which arc in agreement with each other. For example, matching SU(3) onto SU{2) for the 
nucleon mass, one expects D + F = cja + 0(m s ). Furthermore, when the resulting formulae are used to extrapolate to 
the physical point, they are in disagreement with the experimentally measured values of the baryon mass splittings. 
By comparing the predicted NLO contributions to the mass splittings of the octet baryons with their experimentally 
determined values, we demonstrate that a fit that only includes 0(m5), or NLO terms, is doomed to fail to reproduce 
these LECs, and also it is most probable that the SU (3) heavy baryon x?T determination of the mass splittings as 
well as the masses fails to converge. In Table XI, we present the experimental knowledge of the octet baryon masses 
and mass splittings, as well as the predicted NLO contribution to the masses (at the physical point), i.e. we have set 
% = Pm = om = (see Eq. (46)). We have chosen the charge neutral baryon masses, as our lattice calculations 
in which electromagnetism is turned off are most similar to these. To determine the central values and errors of 

the various 5M B , we have used the values of the axial couplings D = 0.715(50) and F = 0.453(50), which were 
taken from Rcf. [24] (we have inflated the uncertainties to be consistent with those in Ref. [68-73]) and we have used 
C = 1-2(2) that is consistent with Rcfs. [9, 69]. We have performed this exercise with heavy baryon xPT both with 
and without (A) explicit decuplet degrees of freedom. 

As can be seen in Table XI, the predicted NLO contributions to the various mass splittings are in all but one case 
larger in magnitude than the actual splitting itself, and generally opposite in sign. To accommodate these NLO mass 
corrections and the experimental results, the LO contributions would have to be larger still and they would not lead to 

(3 /2) (3 /2) (1) (1) 

the expected hierarchy of higher order contributions, SM B , — 5M B <§; SM B , — SM B essential to the convergence 
of the expansion. The explicit inclusion of the decuplet degrees of freedom generally improves the situation but not 
enough to alleviate concerns of convergence. 19 This was discussed in some detail in Rcf. [110]. One exception to this 
picture is the Gell-Mann-Okubo (GMO) formula, where we have defined the mass splittings, 

A GMO = M A + ^M E - ^M N - ?-M s . (54) 

The NLO contributions to the octet baryon masses provide the leading correction to GMO formula, and as can be 
seen in Table XI, the SU (3) HB^PT formula is in remarkable agreement with the experimental measurements. This 



There is a subtle issue that arises with the inclusion of the decuplet states, requiring the inclusion of an extra chiral-singlet parameter 
in the theory, the mass splitting of the decuplet and octet baryons in the chiral limit. Because of this extra parameter, one can 
never completely disentangle the LO and NLO contributions to the baryon masses as the new parameter leads to an arbitrary finite 
renormalization of all the existing LECs. One can therefore only make rigorous statements about the complete mass calculation to a 
given order and not the absolute size of a given order in m q [67, 80, 81]. 
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FIG. 15: The Gell-Mann-Okubo (GMO) mass ratio with 5 G mo = 
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also put the experimental number in the figure, using the charge-neutral baryon masses. 



where Agmo is defined in Eq. (54) and the centroid 
| M A + |Ms + \Mn + jM s . Here we collect the results of this work, as well as those of NPLQCD [14]. We 



can be understood in part from the fact that the quark mass matrix transforms as an 8 1 under the SU (3) flavor 
symmetry. Therefore, a single insertion of the quark mass matrix in the Lagrangian, such as the operators in Eq. (46) 
with coefficients a« and (3m , will automatically satisfy the GMO relation. Violations of the GMO relation must come 
from higher dimensional representations, the first of which is the 27 and appears at NLO in the HB^PT expansion. 
In Ref. [14], the NPLQCD Collaboration provides a more detailed discussion of this topic and presents a detailed 
calculation of the violations of the GMO relation, finding the agreement between the NLO HB^PT formula and the 
lattice calculation holds for heavier quark masses as well. In Figure 15 we display our resulting calculations of the 
deviation of the GMO relation, divided by the centroid octet mass, Mb = \M^ + §Ms + jMjy + jMh, following 
NPLQCD, whose results we find agreement with and are plotted along side ours. 

One must be careful in making this convergence analysis. The leading axial couplings in the Lagrangian, D, F, and 
C are not physical obscrvablcs. However, with sufficiently light up, down and strange quark masses, one would expect 
the fits to give values for these LECs which are fairly stable to the inclusion of higher orders in the chiral expansion. 
These axial couplings provide the leading order contribution to the axial matrix elements that give, gA = D + F, 
9AS = 2D, gsn = 2F and gss = D — F (following the normalization of Ref. [Ill])- To fully test issues of convergence 
for SU(3) HB^PT, one should perform a combined analysis of the octet baryon axial charges and masses, using the 
NNLO expressions for both. This is beyond the scope of this work. We now proceed to perform fits to our octet mass 
results using the HB^PT formula. 



2. Baryon mass extrapolations 



In this section, we perform a large variety of chiral extrapolation analyses to our calculated octet and dccuplet 
baryon masses. We find that, whereas the resulting minimizations are statistically consistent with the lattice results, 
as measured by the x 2 per d.o.f., in general, the resulting chiral extrapolations are in disagreement with the physical 
masses and mass splittings. We perform both the continuum SU (3) analysis as well as the mixed-action analysis, 
using the mass formulae that can be determined from Eq. (46), for example Eq. (48). Before presenting our results, 
we need to note that it is known the strange quark mass on these MILC lattices is too large [112, 113]. As we have 
used only one value for the strange quark mass, we can not control the strange quark mass extrapolation. Attempting 
to correct for this leads to a reduction in all the baryon masses containing strange quarks, which provides for a larger 
disagreement between our extrapolated baryon mass splittings and the physical mass splittings as evidenced from the 
resulting fits presented in Tables XII to XX. 

We begin with a LO two-flavor extrapolation of the octet baryon mass splittings, in which the strange quark is 
integrated out, but we enforce 5(7(3) symmetry. Wc perform the fit to M\ — Mjy, My, — Mjv and Me — Mn using 



32 



M B - M N = otq - a\ 




(a) 

FIG. 16: We display the LO SU(2) fits to the mass splittings, 
mass points plus the 68% confidence bands. The right plot, ( 
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The left plot, (a) is the result of fitting only to the lightest two 
j) is the result of fitting to the lightest three mass points. The 
included in the analysis. 



the following formula 



M A - M N = & s ~ t,®-*™^ , 



Ms - M w = i<5 s + p s - Qc^, + ml , 

M 3 - M N = + p s - f^a u + ^ ml , (55) 

where in terms of the LECs of Eq. (46), we have 

a s = -a M m s , j3 s = -fhim s , a u = — — , fl u = — — . (56) 

-fc>o i>o 

The results of this analysis are collected in Table XII and in Figure 16 we display the resulting chiral extrapolations 
from the analysis using the lightest two and three quark mass ensembles. 

We see that the fits to the lightest two and three quark mass ensembles have a good x 2 per d.o.f., while the inclusion 
of the fourth lightest point, the m030 ensemble, results in a poor fit. Furthermore, these two-flavor extrapolations 



TABLE XII: Results from LO SU{2) extrapolations with St/ (3) symmetry, Eqs. (55) and (56). 



FIT: LO SU{2) 


range 


a s 






0u 


x 2 


d.o.f. 






[MeV] 


[MeV] 


[MeV" 1 ] 


[MeV" 1 ] 






M B - M N 


007-010 


174(13) 


169(15) 


0.00085(22) 


0.00087(25) 


1.1 


2 




007-020 


164(04) 


162(04) 


0.00069(04) 


0.00075(04) 


2.5 


5 




007-030 


158(03) 


150(03) 


0.00061(02) 


0.00061(02) 


42.4 


8 


"Predictions" 


M A - M N 


Me - M N 


M s - Mat M s - M A 


M E - M A 




M B - Ms 


Phys. MeV 


176 


254 


376 


77 


199 




122 


007-010 


166(11) 


270(15) 


384(16) 


104(21) 


218(8) 




114(20) 


007-020 


158(04) 


260(05) 


367(05) 


102(07) 


209(3) 




107(07) 


007-030 


152(03) 


243(04) 


350(04) 


90(05) 


197(2) 




107(05) 
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give mass splittings at the physical pion mass that are reasonably close to the physical octet mass splittings. Next, 
we perform a variety of NLO analysis using both the continuum SU(3) and the mixed-action formulae; 

• M B - M Nl Table XIV 

• m b~ m n Table XV. Dividing by f n removes scale setting ambiguities. 

• 7 K ■ Dividing by Jk also removes scale setting ambiguities, but more importantly fx has significantly 
milder chiral corrections than f^. Therefore, since the chiral corrections to / appear at NNLO in the mass 
splittings, beyond the order of this fit, using fx instead of /„■ provides more stable fits, as seen in Table XV. 

• M B - M N and M N , Table XVI. 

• M N , M A , Afs and Ms, Table XVII. 

• Mb — M/v and M/v at NNLO. We do not find a stable minimization for this fit, as there are 18 unknown LECs 
that must be determined. 

• Mb = a B + afm^, the straight line analysis as for the nucleon, Table XIII. 

In Figure 17 we display some of the resulting extrapolations. In all cases, the (gray) boxes denote points that were 
not included in the minimization. The (colored) circles with error bars are points that were included, and the dashed 
boxes with error bars arc the resulting predictions at the physical pion mass, slightly displaced horizontally for clarity. 
The (colored) stars are the physical masses (mass splittings) that are never included in the minimization analysis. 
We perform a similar analysis for the dccuplct baryon masses, with the same caveat mentioned for the delta mass 
in Sec. IV B 3. Aside from the O - , the other decuplet states have similar volume issues as the delta discussed in 
Sec. IV B 3. The resulting fits are collected in Tables XVIII to XX. 

There are a few points worth mentioning. First, as measured by the x 2 per d.o.f., most of the fits presented in 
Tables XII to XX provide a good description of our lattice results, and are in disagreement with the physical baryon 
masses and mass splittings. One should not place much emphasis on the disagreement with the physical masses. 
First, as mentioned above, the strange quark mass is known to be too large [112, 113]. Furthermore, the issues of 
convergence of SU(3) HB^PT may lead to large NNLO corrections, or worse, indicate a lack of convergence. This 
problem has been discussed in some detail in Ref. [76, 114]. Support of this statement is found in the LO analysis 
presented in Table XII (and displayed in Fig. 16). The resulting LECs and predicted baryon masses (mass splittings) 
from the various fits are all in reasonable agreement with each other, despite providing axial couplings, D, F, C and 
H that are in stark disagreement with their known phenomenological values. We should mention that the inclusion 
of the NNLO terms of 0{m\), may provide stability to the fits, such that the values of these couplings are in closer 
agreement with phenomenology. However, performing the NNLO analysis with D, F, C and H fixed, we also do not 
find a stable minimization. 20 Examining Tables XVIII through XX, it is clear that the decuplet extrapolation is in far 
worse condition than the octet masses. Understanding the decuplet masses will require a multiple volume study, which 
is beyond the scope of this work. Lastly, we note that for the straight line in analysis, using Mb = a B + afm„, 
we find (see Table XIII) that the resulting values of the parameters af are approximately simple fractions: ~ 1, 
cti ~ 2/3, af ~ 1/2 and af ~ 1/3. Determining whether this is a phenomenon of QCD or perhaps a combined 
finite-volume-lattice-spacing artifact will require further investigation with multiple lattice spacings and volumes. 



TABLE XIII: Results from straight-line in m-n fit of the octet masses. A noteworthy feature is the coefficient of the m-n term 
for each of the masses. 



M B = a B + afm^ 


M N 


M A M s 


Me 


Mjv 


M A 


M E 


M s 






a,f[GeV] 








af 




007-030 


0.82(2) 


1.05(1) 1.19(1) 


1.33(1) 


0.94(4) 


0.65(3) 


0.45(3) 


0.29(2) 


007-040 


0.80(1) 


1.03(1) 1-17(1) 


1.32(1) 


0.99(3) 


0.69(2) 


0.50(2) 


0.32(2) 


007-050 


0.79(1) 


1.03(1) 1.16(1) 


1.31(1) 


1.01(2) 


0.70(2) 


0.52(2) 


0.33(1) 



Holding H fixed may not be ideal regardless, as in a recent investigation of H in SU(2) xPT (matching provides H = 9aa)j it was 
found that H has a large expected quark mass dependence [115]. 
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FIG. 17: We display various chiral extrapolations of the octet baryon masses and mass splittings. In all fits, the (gray) squares 
are points not included in the analysis, the (colored) open circles with error bars are the points that are included and the 
(colored) stars are the physical masses/splittings, which are never included in the minimization. The dashed (colored) squares 
with error bars are the resulting predictions, slightly displaced horizontally for clarity. Figs, (a) and (b) are from a combined 
fit of Mb — Mn and Mn, fit to the lightest three quark mass ensembles, using the NLO mixed-action extrapolation formulae. 
Fig. (c) is the result of the NLO mixed-action fit to the octet baryon masses. Fig. (d) is the result of the straight line analysis. 
As can be seen, the straight line fit reproduces the lattice results not included in the minimization but conflicts dramatically 
with experiment. The error bars/bands represent the 68% confidence interval from the statistical uncertainty. 



V. CONCLUSIONS 



In this work, we have presented a detailed study of the light hadron spectrum in the mixed action framework 
that consists of computing domain-wall valence fcrmion propagators on the background of the Asqtad improved, 
rooted, staggered MILC sea configurations. These results assume the validity of the rooting procedure used in the 
staggered sea sector. With the domain wall pion masses tuned to within a few percent of the staggered Goldstone pion 
masses, we found that the other meson masses had no systematic trend compared to the masses determined by MILC. 
However, as discussed in detail in Sec. IIIB, there are additional complications which must be addressed to compare 
the vector meson masses computed in this work to those of MILC. In contrast, we found that the baryon masses 
were systematically lighter (and therefore in better agreement with the experimental values) than those computed on 
the coarse MILC lattices, suggesting the mixed action has smaller lattice spacing corrections. This is supported in 
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TABLE XIV: Results from NLO bootstrap ^-extrapolations of the octet mass splittings, using mixed action (MA) and SU(3) 
heavy baryon xPT- 



FIT- 1SJT n 


TcLIlgG 




O-M 


PM 


r< 




r 


X 


j n f 

Q.O.I. 








[GeV- 1 ] 


[GeV- 1 ] 












Mb - M N 


007-020 


MA 


-0.49(04) 


-0.43(03) 


0.38(3) 


0.07(07) 


0.03(3) 


10.1 


4 




007-020 


SU{3) 


-0.87(38) 


-0.40(18) 


0.36(3) 


0.23(19) 


0.27(9) 


3.1 


4 




007-030 


MA 


-0.48(03) 


-0.43(03) 


0.38(2) 


0.06(07) 


0.02(3) 


11.4 


7 




007-030 


SU(3) 


-0.77(32) 


-0.32(19) 


0.37(2) 


0.15(19) 


0.27(8) 


4.8 


7 




007-040 


MA 


-0.48(04) 


-0.44(04) 


0.38(2) 


0.10(08) 


0.03(3) 


11.6 


10 




007-040 


SU(S) 


-0.86(36) 


-0.38(20) 


0.37(2) 


0.21(20) 


0.28(8) 


8.8 


10 


"Predictions" 


M A - M N 


M s - Mjv 


M E 


-M N 


M E - M A 


Ms - M A 




M= - M E 


Phys. MeV 




176 


254 




376 


77 


199 




122 


007-020: MA 




115(3) 


182(4) 


262(4) 


68(5) 


148(3) 




80(4) 


007-020: SU(3) 




131(13) 


195(17) 


286(27) 


64(4) 


154(14) 




90(11) 


007-030: MA 




115(3) 


182(4) 


262(4) 


68(4) 


147(3) 




79(4) 


007-030: SU(3) 




127(9) 


193(14) 


276(19) 


66(7) 


149(11) 




83(5) 



the case of the nucleon mass by comparing to the (preliminary) super-fine MILC results displayed in Fig. 12(a). We 
also find our nucleon mass results are consistent with those calculated by the RBC/UKQCD group using a domain 
wall valence on domain wall sea action, as can be inferred from Fig. 12(a) and (c), as well as with the results of the 
QCDSF/UKQCD group using two-flavors of 0(a) improved Wilson fermions, inferred from Fig. 12(6). 

In addition to the spectroscopy calculations, we have performed a detailed chiral extrapolation analysis of the 
octet and decuplet baryon masses, using both the continuum SU(3) heavy baryon xPT as well as its mixed action 
generalization. We have performed extensive three flavor chiral extrapolations, the results of which are collected in 
Tables XIV-XX. In most cases, the extrapolation analysis is in good agreement with our calculational results, as 
measured by the x 2 per degree of freedom. However, in all cases, the resulting values of the axial couplings, D and F 
are different from those determined phcnomcnologically. This is suggestive that the three flavor heavy baryon xPT 
is not converging for these masses. Furthermore, the analysis presented in Table XI is suggestive that the theory is 
not converging even at the physical point, due to the large kaon mass, except for special observables like the Gell- 
Mann-Okubo relation. To form a conclusive analysis of the convergence SU(3) heavy baryon xPT, one would need to 
perform a chiral extrapolation of the octet masses (or mass splittings) simultaneously with the hyperon axial charges, 
with enough lattice-data points to perform the NNLO analysis. This is beyond the scope of this work. Alternatively, 
one can perform SU(2) chiral extrapolations of the masses, as advocated in Rcf. [116], which will have much better 
convergence properties. We have performed the leading order SU (2) extrapolation of the hyperon mass splittings, 
while enforcing SU (3) symmetry, the results of which are collected in Table XII. Wc found good agreement with both 
the our numerical results as well as the physical mass splittings. 

We have also performed the continuum SU (2) chiral extrapolation of the nucleon and delta masses using the O(m^) 
mass formulae. These extrapolations required the input of the nucleon and delta axial couplings, gA, <?a a and <?atv- 
Even after fixing these parameters, the mixed action formula has too many unknown parameters to be determined 
from our lattice results for either the nucleon or delta, and so we did not perform this mixed action analysis. In 
the case of the delta mass, the extrapolation does not agree with the known pole position of the A(1232), even 
with the inclusion of the predicted large volume corrections from Ref. [103]. Furthermore, the delta mass on the 
one larger volume (L ~ 3.5 fm) ensemble we use is heavier than on the smaller volume (L ~ 2.5 fm), yielding a 
stronger discrepancy with the physical pole mass. For the lightest two mass points, on one of which we have the two 
volumes, the delta would decay in infinite volume but due to the restricted values of allowed momenta, this decay is 
kincmatically forbidden. From our analysis, it is clear that an understanding of the delta mass as well as its chiral 
extrapolation will require a multi-volume study where the finite volume systematics can be carefully explored. Wc 
expect similar issues for the S* and S* decuplet baryon resonances. 

We found that the 0(m*) SU(2) chiral extrapolation of the nucleon mass is in statistical agreement with the lattice 
data as well as the physical nucleon mass. An examination of resulting contributions from the first three orders in the 
expansion, however, revealed that only for the lightest one or two of six mass points does one trust the resulting fits to 
the nucleon mass in the strict sense of an order by order convergence. The most surprising result of our analysis is that 
a linear fit in the pion mass, = Oq -\-a^m v , Eq. (22), provides a remarkable agreement with both the lattice data 
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TABLE XV: Results from NLO extrapolations of the octet mass splittings in f w and /k units, using mixed action (MA) and 
SU{3) heavy baryon xPT. 



PTT- NT O 

J? J. _L . IN 1 j V .7 






l bev J 


flu r 

I'M 

l bev J 


r< 


Y) 


p 


X 


H n f 

Ll.U.l. 


M B -M N 
h 


007-020 


MA 


-1.7(4) 


-0.8(2) 


0.00(41) 


0.50(05) 


0.33(09) 


2.2 


4 




007-020 


SU(3) 


-0.6(2) 


-0.5(3) 


0.27(08) 


0.09(80) 


0.04(26) 


3.1 


4 




007-030 


MA 


-0.6(8) 


-0.5(2) 


0.25(10) 


0.08(51) 


0.10(70) 


5.1 


7 




007-030 


SU(3) 


-0.6(2) 


-0.6(2) 


0.25(06) 


0.20(24) 


0.00(16) 


4.8 


7 




007-040 


MA 


-0.6(3) 


-0.5(4) 


0.23(08) 


0.1(1.1) 


0.02(35) 


10.6 


10 




007-040 


SU(3) 


-0.6(2) 


-0.6(2) 


0.21(07) 


0.32(15) 


0.00(10) 


8.8 


10 



"Predictions" 


Ma - M N 


Ms - M N 


M s - Mjv 


M E - M A 


M H - M A 


Ms - Me 


Phys. MeV 


176 


254 


376 


77 


199 


122 


007-020: MA 


220(43) 


233(09) 


443(55) 


13(48) 


223(13) 


210(58) 


007-020: 5*17(3) 


134(8) 


209(16) 


306(21) 


74(14) 


171(16) 


97(09) 


007-030: MA 


139(61) 


210(14) 


313(93) 


71(50) 


174(33) 


103(81) 


007-030: SU(3) 


133(8) 


210(09) 


306(14) 


78(10) 


173(09) 


95(09) 



range 




CtM 

[GeW 1 ] 


(3m 
[GeW 1 ] 


C 


D 


F 


x 2 


d.o. 


007-020 


MA 


-1.1(2) 


-0.8(2) 


0.14(19) 


0.57(13) 


0.24(06) 


4.0 


4 


007-020 


SU(3) 


-0.8(1) 


-0.8(2) 


0.19(14) 


0.53(13) 


0.18(04) 


5.4 


4 


007-030 


MA 


-0.9(2) 


-0.8(2) 


0.24(10) 


0.48(14) 


0.18(05) 


7.4 


7 


007-030 


SU(3) 


-0.7(1) 


-0.7(2) 


0.26(10) 


0.45(13) 


0.14(04) 


8.0 


7 


007-040 


MA 


-0.8(1) 


-0.9(2) 


0.19(12) 


0.52(14) 


0.16(04) 


12.4 


10 


007-040 


517(3) 


-0.7(1) 


-0.8(2) 


0.20(11) 


0.48(12) 


0.12(03) 


12.9 


10 



"Predictions" Ma - M N Me - M N Ms - M N Me - Ma Mh - Ma Mh - Me 

Phys. MeV 176 254 376 77 199 122 

007-020: MA 168(20) 234(25) 373(32) 66(33) 205(19) 139(29) 

007-020: 5(7(3) 136(07) 221(14) 321(17) 85(13) 185(13) 100(07) 

007-030: MA 114(15) 223(23) 338(26) 74(27) 189(17) 115(22) 

007-030: 5(7(3) 131(06) 211(12) 307(14) 80(12) 175(11) 95(06) 



as well as the physical nucleoli mass. To explore our uncontrolled systematics, the lattice spacing and finite volume 
artifacts, we compared our results to the recent nucleon mass calculations of the MILC [96], QCDSF/UKQCD [93], 
RBC/UKQCD [97] and ETM [94] Collaborations. We found that Eq. (22) provides a statistically good description 
of the nucleon mass lattice data of all these groups, although the heavy baryon extrapolation was noticeably better 
in the case of the ETM calculation. This straight line in m n analysis is not based upon any understanding of low 
energy QCD we currently have; it is a lattice phenomenological form. At this point, we are unable to determine if 
this phenomenon is a conspiracy of QCD or whether it arises from lattice artifacts. To resolve this issue, more lattice 
calculations are needed at lighter pion masses, as well as multiple lattice spacings and volumes. The RBC/UKQCD 
Collaborations, which use 2 + 1 flavor domain- wall valence and sea fermions with a ~ 0.114 fm, have nucleon mass 
results that are consistent with ours. The MILC Collaboration has performed high statistics calculations with three 
lattice spacings ranging from the coarse lattice, a ~ 0.124 fm, to the super-fine lattices with a ~ 0.06 fm. While the 
coarse MILC nucleon mass results are systematically higher than ours, the super-fine nucleon mass results, with a 
pion mass as light as m v ~ 220 MeV, lie within the statistical errors on the straight line drawn through our mass 
results, as can be seen in Fig. 12. As is clear from this analysis, a reliable chiral extrapolation of the nucleon mass 
will be much more challenging and demanding than perhaps previously thought. The nucleon mass displays too little 
structure to determine all the LECs from an extrapolation of the nucleon mass alone. To make progress, one will 
need to perform a global extrapolation analysis with a sufficient set of observables designed to strongly constrain the 
values of the LECs, notedly the axial couplings, gA, g&N and <?aa- We leave this to future work. 
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TABLE XVI: Results from NLO bootstrap ^-extrapolations of the octet mass splittings combined with the nucleon mass, using 
mixed action (MA) and SU(3) heavy baryon x?T- 



FIT: NLO 


ranf 


;e 




Mo 


om 




Pm 


C 


D 


F 


x 2 


d.o.f. 










[GeV] 


[GeV" 1 ] 


[GeV" 1 ] 


[GeV" 1 ] 












M B - M N 


007- 


-020 


MA 


0.807(50) 


-0.23(5) 


-0.49(3) 


-0.43(2) 


0.37(3) 


0.07(6) 


0.03(3) 


10.9 


5 


M N 


007- 


-020 


SU{3) 


0.895(43) 


-0.25(6) 


-0.80(26) 


-0.37(11) 


0.36(2) 


0.21(13) 


0.26(7) 


9.2 


5 




007- 


-030 


MA 


0.848(38) 


-0.19(3) 


-0.49(3) 


-0.43(2) 


0.38(2) 


0.09(6) 


0.03(3) 


16.8 


9 




007- 


-030 


517(3) 


0.952(66) 


-0.21(7) 


-0.94(38) 


-0.42(22) 


0.36(2) 


0.27(20) 


0.29(8) 


13.1 


9 




007- 


-040 


MA 


0.864(30) 


-0.18(3) 


-0.50(4) 


-0.45(3) 


0.38(2) 


0.14(5) 


0.04(3) 


17.0 


13 




007- 


-040 


517(3) 


0.946(55) 


-0.22(6) 


-0.92(36) 


-0.41(20) 


0.36(2) 


0.26(20) 


0.29(8) 


13.3 


13 



"Predictions" 


M A - Mat 


M E - Mat 


M s - M N 


Me - Ma 


Ms - M A 


Ms - Mr, 


Mat 


M A 


Me 


Ms 


Phys. MeV 


176 


254 


376 


77 


199 


122 


938 


1116 


1193 


1314 


007-020: MA 


115(02) 


182(03) 


262(04) 


67(4) 


147(2) 


81(3) 


952(24) 


1066(24) 


1133(23) 


1214(24) 


007-020: SZ7(3) 


129(10) 


185(07) 


271(14) 


60(2) 


147(7) 


86(7) 


980(20) 


1105(23) 


1165(23) 


1251(29) 


007-030: MA 


115(02) 


182(03) 


262(04) 


67(3) 


148(2) 


80(3) 


968(20) 


1083(19) 


1150(18) 


1231(19) 


007-030: 5(7(3) 


125(07) 


196(18) 


281(23) 


68(7) 


153(13) 


85(6) 


1004(25) 


1133(33) 


1200(40) 


1285(45) 



TABLE XVII: Results from NLO bootstrap ^-extrapolations of the octet baryon masses, using mixed action (MA) and 5(7 (3) 
heavy baryon xPT- 



FIT: NLO 


ranj; 


;e 




Mo 


cm 


CtM 


Pm 


C 


D 


F 




d.o.f. 










[GeV] 


[GeV- 1 ] 


[GeV" 1 ] 


[GeV- 1 ] 












Mat, Ma, 


007- 


-020 


MA 


1.087(51) 


-0.03(5) 


-0.72(8) 


-0.62(4) 


0.15(9) 


0.33(4) 


0.14(3) 


6.0 


5 


Me, Ms 


007- 


-020 


5(7(3) 


1.014(32) 


-0.07(4) 


-0.77(10) 


-0.56(5) 


0.18(9) 


0.30(6) 


0.19(4) 


5.5 


5 




007- 


-030 


MA 


1.149(57) 


0.01(4) 


-0.79(11) 


-0.67(7) 


0.12(9) 


0.38(6) 


0.16(3) 


14.4 


9 




007- 


-030 


5(7(3) 


1.091(66) 


-0.04(3) 


-0.99(28) 


-0.73(19) 


0.1(1) 


0.44(14) 


0.24(7) 


11.9 


9 




007- 


-040 


MA 


1.147(52) 


0.01(3) 


-0.78(10) 


-0.68(6) 


0.13(9) 


0.39(6) 


0.16(3) 


14.9 


13 




007- 


-040 


5(7(3) 


1.090(61) 


-0.04(3) 


-0.99(26) 


-0.73(18) 


0.1(1) 


0.45(13) 


0.25(6) 


12.5 


13 



"Predictions" 


M A - Mat 


Me - M N 


M s - M N 


Me - M A 


Ms - M A 


Ms - Me 


M N 


M A 


Me 


Ms 


Phys. MeV 


176 


254 


376 


77 


199 


122 


938 


1116 


1193 


1314 


007-020: MA 


122(4) 


183(3) 


277(6) 


61(6) 


155(2) 


94(7) 


1082(23) 


1204(26) 


1266(23) 


1360(27) 


007-020: 5(7(3) 


117(3) 


181(2) 


265(4) 


65(3) 


148(2) 


83(3) 


1028(17) 


1145(16) 


1240(23) 


1293(17) 


007-030: MA 


125(5) 


184(3) 


283(8) 


59(6) 


158(4) 


99(8) 


1112(27) 


1237(32) 


1296(28) 


1395(34) 


007-030: 5(7(3) 


121(8) 


189(9) 


277(16) 


68(3) 


156(8) 


88(7) 


1051(18) 


1172(22) 


1237(21) 


1328(29) 
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TABLE XVIII: Results from NLO bootstrap X" ex trapolations of the decuplet mass splittings using mixed action (MA) and 
SU(3) heavy baryon xPT. 



PTT^- TvTT O 
rll. 1\LU 


range 




7a/ 


L 


u 

±1 


■v 2 


a. o.i. 








[GeV] 


[GeV" 1 ] 


[GeV" 1 ] 






M T - Ma 


007-020 


MA 


1.2(2) 


0.00(14) 


1.19(32) 


8.4 


6 




007-020 


8U(3) 


0.5(6) 


0.83(14) 


0.00(93) 


3.4 


6 




007-030 


MA 


0.48(01) 


0.00(18) 


0.75(35) 


11.7 


9 




007-030 


8U(3) 


0.5(5) 


0.65(15) 


0.00(82) 


8.1 


9 




007-040 


MA 


0.95(13) 


0.00(15) 


0.99(21) 


13.2 


12 




007-040 


SU(3) 


0.5(5) 


0.70(11) 


0.00(86) 


8.3 


12 


"Predictions" 






M E . - M A 




M s . - Ma 




M n - - M A 


Phys. MeV 






152 




300 




440 


007-020: MA 






110(20) 




207(32) 




292(38) 


007-020: SU(3) 






73(03) 




142(05) 




207(08) 


007-030: MA 






87(13) 




169(21) 




246(25) 


007-030: 5(7(3) 






71(02) 




140(05) 




206(08) 



TABLE XIX: Results from NLO bootstrap x -ex trapolations of the decuplet mass splittings combined with the delta mass, 
using mixed action (MA) and 5(7(3) heavy baryon x?T. 



FIT: NLO 


ranj 


;e 




Mr,o 


cm 




C 


H 


X 2 d.o.f. 










[GeV] 


[GeV" 1 ] 


[GeV" 1 ] 










Mr - Ma 


007- 


-020 


MA 


1.74(14) 


-0.01(08) 


1.2(4) 
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1.22(31) 


11.4 


7 


Ma 


007- 


-020 


5(7(3) 


1.69(12) 


0.18(11) 


0.50(02) 


0.71(16) 


0.00(1.02) 


10.3 


7 




007- 


-030 


MA 


1.55(10) 


-0.05(05) 


0.8(3) 


0.00(16) 


0.80(33) 


15.5 


11 




007- 


-030 


577(3) 


1.54(08) 


0.03(07) 


0.48(01) 


0.52(18) 


0.00(80) 


14.7 


11 




007- 


-040 


MA 


1.66(08) 


0.00(03) 


0.9(2) 


0.00(30) 


0.99(20) 


18.8 


15 




007- 


-040 


5(7(3) 


1.62(07) 


0.10(05) 


0.49(01) 


0.57(12) 


0.00(1.04) 


19.2 


15 


"Predictions" 




M E 


. - M A 


M=« - M A 


M Q - - Ma 


M A 


Me* 


M s . 


M n - 




Phys. MeV 
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1232 


1384 


1532 
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007-020: MA 
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210(32) 


295(38) 


1627(75) 


1739(92) 
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007-030: MA 




88(13) 
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1533(55) 


1622(66) 
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1782(77) 




007-030: 577(3) 




72(02) 


142(05) 


211(8) 


1520(40) 


1592(40) 


1662(39) 


1731(37) 
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